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ABSTRACT 

We present an updated and revised analysis of the relationship between the H/3 broad-line region (BLR) ra- 
dius and the luminosity of the active galactic nucleus (AGN). Specifically, we have carried out two-dimensional 
surface brightness decompositions of the host galaxies of 9 new AGNs imaged with the Hubble Space Tele- 
scope Wide Field Camera 3. The surface brightness decompositions allow us to create "AGN-free" images 
of the galaxies, from which we measure the starlight contribution to the optical luminosity measured through 
the ground-based spectroscopic aperture. We also incorporate 20 new reverberation-mapping measurements of 
the H/3 time lag, which is assumed to yield the average H/3 BLR radius. The final sample includes 41 AGNs 
covering four orders of magnitude in luminosity. The additions and updates incorporated here primarily affect 
the low-luminosity end of the Rblr-L relationship. The best fit to the relationship using a Bayesian analysis 



finds a slope of a = 0.533; 



0.033' 



consistent with previous work and with simple photoionization arguments. 



Only two AGNs appear to be outliers from the relationship, but both of them have monitoring light curves that 
raise doubt regarding the accuracy of their reported time lags. The scatter around the relationship is found to 
be 0. 19 ± 0.02 dex, but would be decreased to 0. 13 dex by the removal of these two suspect measurements. A 
large fraction of the remaining scatter in the relationship is likely due to the inaccurate distances to the AGN 
host galaxies. Our results help support the possibility that the Rblr-L relationship could potentially be used 
to turn the BLRs of AGNs into standardizable candles. This would allow the cosmological expansion of the 
Universe to be probed by a separate population of objects, and over a larger range of redshifts. 
Subject headings: galaxies: active — galaxies: nuclei — galaxies: photometry — galaxies: Seyfert 
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1. INTRODUCTION 
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The ability to determine black hole masses is a crucial 
step toward understanding the link between galaxies and their 
black holes, as well as the details of the black hole environ- 
ment. To date, dynamical methods have resulted in measure- 
ments of some 50 black hole masses. However, these meth- 
ods require that the gravitational influence of the black hole 
on the stars or gas be spatially resolved, effectively limiting 
the reach of current dynamical methods to galaxies no further 
than ~ 150Mpc for even the most massive black holes (see 
iGultekinet alj2009l) . 

Active galactic nuclei (AGNs), in contrast, are some of the 
most luminous objects in the Universe, and are thus capable of 
providing us with the leverage needed to probe the growth and 
evolution of black holes at any significant cosmological dis- 
tance. However, AGNs are also so rare that even the nearest 
are generally too distant for current instruments to spatially 
resolve the radius of influence of the black hole and provide 
a local calibration for their masses. Instead, the most suc- 
cessful technique for meas uring black hole masses in AGNs 
is reverberation mapping (Blandford & McKee 1982). Re- 
verberation mapping requires high-quality spectrophotomet- 
ric monitoring of an AGN over an extended period of time. 
The line-emitting regions that give rise to the characteristic 
AGN spectral signatures are photoionized by the hot accretion 
disk around the black hole. The continuum flux (which arises 
from the accretion disk or very close to it) varies with time, 
and these variations are echoed later by changes in the flux of 
the broad emission lines. The delay time between the contin- 
uum variations and the broad-line variations can be measured 
by cross correlation of the light curves and gives the light- 
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travel time across the broad-line region (BLR), or the radius of 
the BLR when multiplied by the speed of light. In effect, re- 
verberation mapping substitutes high temporal resolution for 
high spatial resolution, allowing us to probe regions of gas 
that are only ~ 0.01 pc in extent (comparable to t he inner re- 
gion o f the Oort cloud in our own Solar System; iBrown et al.l 
2004) in the centers of arbitrarily distant galaxies. Combin- 
ing the BLR radius with the mean velocity of the BLR gas, as 
measured from the Doppler broadening of the emission lines, 
and assumptions or indirect estimates of the virial coefficient 
gives a direct constraint on the black hole mass via the virial 
theorem. 

The validity of reverberation masses has been upheld by 
several independent lines of evidence. A subset of objects 
in the reverberation sample have measurements for several 
different emission lines throughout the ultraviolet and opti- 
cal portions of their spectra, and the multiple emission lines 
show a virial behavior (e.g.. iPetersqn & Wandelll999[ 120001 
lKollatschny||2003l: iBentz et alj|2010h . Two AGNs in the cur- 
rent reverberation sample — NGC 3227 and NGC 415 1 — are 
sufficiently close that dynamical modeling has successfully 
determined their black hole masses, and both the stella r dy- 
namical masses (Davies et al. 2006; Onken et al J 120071) and 
gas dynamical masses dHicks & Malkan 2008) agree with the 
reverberation-based masses within the uncertainties. Further- 
more, a fully general Bayesian modeling code has recently 
been developed to analyze reverberation-mapping datasets 
and place limits on th e black hole mass a nd the BLR ge- 
ometry and dynamics dPancoast et al.ll201rl) . When a pplied 
to the reverberatio n mapping data for A rp 151 dBrewer et al.l 
120111) and Mrk50 dPancoast et al.ll2012l) . the method recovers 
a black hole mass that is essentially the same value as that 
determined from the reverberation method outli ned above, 
for standard assumptions of the virial coefficient (Bentz et al. 
2009b; iBarth et alJEoTTl) . 

Reverberation ma pping has yielded black hole masses fo r 
~ 50 AGNs thus far ([Peterson et al.l2004tlBentz et al.l2 009b). 
The BLR radius-luminosity correlation (/?blr oc L a ) derived 
from this reverberation sample is the basis for all secondary 
techniques used to es t imate black hole masses in distant 
AGNs (e.g.,lLaorfll998l:IWandel et al.l ll999l:lMcLur e & Jarvisl 
2002; Vestergaard & Peterson 2006). The power of the ^blr- 
L relationship comes from the simplicity of using it to quickly 
estimate Mbh for large samples of objects, even at high red- 
shift, with only a single spectrum per object. 

This simplicity has led to the Rblr-L relationship be- 
ing heavily utilized in the literature. A small sampling 
of studies that have utilized the R^lr-L relationship in 
the last few years inclu des investigations of Mbh in the 
most distant quasars (e.g., IWillott et al.ll20Tot iMortlock et all 
l20Tlh . black hole mass functions and Eddin gton ratio 
distributions through cosmic history (e.g.. iGreene & Hoi 
2007bl: l\^stergaard etal.l 120081: IWs tergaard & Osmerl 120091 
Kelly et alJl2009l: ISchulze & WisotzkilfeoiOh , cos mic evolu- 
tion of black holes and the i r host galaxies (e.g., IWoo et all 
120081: iMerloni et all l20Tot iBennert et al.l 12010). the na- 
ture o f narrow-line Seyfert 1 galaxies (e.g., iMathur et al.l 
20121: Papadaki s et al] 1201 Oh . duty cycles of quasars (e.g., 
Shank ar et al.l l2009|). accretion pr o pertie s of v arious types 
of AG Ns (e.g.. IWang etaf] 120091: lCaol I20T0I: ITrump etal.1 
1201 ll) . studies of relativistic jets and the jet-disk con- 
nection (e.g.. ISambruna et ail l2006t iTavecchio et al .1 120071: 
Ghisellini & Tavecch io 2010), studies of black holes in dwarf 
and low-mass galaxies (e.g jGreene & Holl2007al:lbong et al] 



120071: iThornton et al]|2008l) . and studies of optical transients 
(e.g.. lDrakeetal]|2011l) . 

Because of the utility of the Rblr-L relationship, much 
work has gone into removing biases and noise from the rever- 
beration database. Previous determinations of the R^lr-L re- 
lationship used luminosity measurements f rom ground-based 
spect ra and found the slope to be a « 0.7 (Kaspi et al] |2000l 
2005). To achieve the low level of uncertainties in the flux 
calibration necessary for reverberation mapping, a large (e.g., 
5" x 7") spectroscopic aperture is typically employed. There- 
fore, for all the nearby objects with reverberation masses, 
a substantial fraction of the observed luminosity is actually 
the result of the host-galaxy starlight and not the AGN itself. 
The entire low-luminosity end of the R^lr-L relationship was 
o verestimated in L, in effect, artificially steepening the slope. 

Bent z et al.l (|2006a) analyzed HST ACS images of the near- 
est reverberation-mapped AGNs and their host galaxies taken 
through the F550M medium-band V filter. The flux contri- 
bution of starlight through the ground-based spectroscopic 
monitoring aperture was measured from each image, and the 
reverberation-mapping luminosities were corrected accord- 
ingly. The resultant R^lr-L relationship was found, as ex- 
pected, to have a much flatter slope (a = 0.52 ±0.04 compared 
to a = 0.67 ± 0.05), consistent with simple photoionization 
expectations. Consequently, all of the remaining objects in 
the reverberation-mapped sample were imaged with ACS in a 
similar manner with the intent of properly accounting for the 
starlight in each object, even when that contribution was as- 
sumed to be small. While the s lope of the relation ship did not 
change much (a = 0.52 ± 0.06: lBentz et al J200 9a). the scatter 
in the relationship was reduced from ~ 40% to ~ 35%. The 
scarcity of measurements anchoring the low-luminosity end, 
in particular, then became apparent. 

In the meantime, much effort has gone into replacing noisy 
and poorly sampled reverberation datasets and increasing the 
overall range of BLR radii probed. The last several years 
in particular have seen a huge amount of effort invested in 
reverberation-mapping experiments that preferentially target 
AGNs with relatively low luminosities. The Lick AGN Mon- 
itoring Project (LAMP) campaign targeted low-luminosity 
AGNs to more fully populate the low-luminosity end of the 
Rblr-L relationship an d succeeded in me asuring H/3 BLR 
radii for 8 new AGNs (IBentz et al]|2009bl) . Multiple recent 
campaigns at MDM Observatory have mainly focused on re- 
placing poorly sampled or noisy reverberation datasets with 
high signal -to-noise ratio (S/N), high temporal c adence spec- 
troscopy (Denney et al. 2 10t iGrier et al.ll20 1 2al) to allow bet- 
ter constraints on the BLR structure and kinematics. 

Given the number of updates, improvements, and additions 
to the reverberation database, we undertook a full recalibra- 
tion of the R^lr-L relationship in an effort to provide a more 
accurately calibrated relationship for the community to use 
when estimating black hole masses in AGNs. In particular, 
our new calibration is more accurate at the low-luminosity 
end where L„ and sub-L* galaxies tend to reside. 

We assume a standard ACDM cosmology of Ho = 
72kms _1 Mpc" 1 , J7m = 0.3, SIa = 0.7 throughout this work. 

2. NEW H/3 BLR MEASUREMENTS 

Recent reverberation-mapping campaigns have focused 
mainly on the low-luminosity end of the R^lr-L relationship 
and provide several new H/3 BLR measurements to the rever- 
beration sample. The measurements come in three separate 
flavors: (1) replacement measurements for targets of previous 
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reverberation campaigns for which the light curves were noi sy 
or undersampled and led to poor or biased determinations of 
the H/3 BLR radius, (2) additional measurements for targets of 
previous reverberation campaigns that already have accurate 
H/3 time-lag measurements, and (3) H/3 BLR measurements 
for new objects that have not been previously examined with 
reverberation mapping. We provide a brief summary of each 
of the monitoring programs with the new results that we incor- 
porate here. The interested reader should refer to the original 
manuscripts reporting the H/3 BLR measurements for more 
details. 

MDM 2005— 3C 390.3 was the subject of a 2005 monitor- 
ing campaign at MDM Observatory that resulted in an addi- 
tional H/3 radius and luminosity measurement for this object 
dDietrichet al.ll2012l) . 

MDM 2007— iDennev et all d2010l) describe the results of a 
2007 monitoring program at MDM and other observatories 
that had a goal of obtaining high-quality, densely sampled 
light curves to search for velocity-resolved time lags in the 
emission lines. H/3 BLR measurements were derived for six 
AGNs through this program, three of which were replace- 
ments for poor-quality datasets, and two of which were ad- 
ditional measurements for AGNs with other reliable measure- 
ments. The final object, Mrk290, was new. Mrk290 was also 
included in the 2008 LAMP sample of AGNs (see below) but 
did not exhibit strong variations during that campaign. The 
inclusion of Mrk 290 in the LAMP sample led to it being in- 
cluded in the HST Cycle 17 imaging campaign that we de- 
scribe below, from which we are able to derive the starlight 
correction to the luminosity. 

LAMP 2008— The 2008 Lick AGN Monitoring Project 
(LAMP) targeted AGNs with estimated black hole masses 
in the range 1O 6 -1O 7 M . Measurements of the H/3 BLR 
radius were determined for eight new objects, and an addi- 
tional measurement of the H/3 BLR radi us in the well-studie d 
AGN NGC5548 was also determined dBentz et alj |2009b). 
NGC 5548 was the only galaxy in the LAMP sample with the 
appropriate HST imaging to allow a host-galaxy starlight cor- 
rection. In the next sections we detail the HST Cycle 17 imag- 
ing program through which we obtained the necessary im- 
ages for the remainder of the LAMP sample, the host-galaxy 
surface brightness modeling of those images, and the de- 
rived starlight corrections to the ground-based spectroscopic 
monochromatic luminosities at 5100 A. 

MDM 2010 — Additional measurements of H/3 radii were 
determined for four AGNs in the reverberation sample dur- 
ing a 2010-201 1 campaign at MDM and other observatories 
dGrieret al.N2012at iPeterson et ail 120 13b . Two other AGNs, 
Mrk6 and Mrk 1501, were new targets and reliable H/3 radii 
were determined for them. Unfortunately, there is no suitable 
HST imaging from which to measure the starlight correction 
to the spectroscopic luminosity. We are therefore unable to 
include them in this analysis of the Rrlr-L relationship. 

The addition of nine new AGNs to the reverberation sam- 
ple along with 1 1 replacement or additional datasets for pre- 
viously monitored AGNs allows us to revisit the calibration of 
the Rblr-L relationship, and in particular to examine the form 
of the relationship at the lower luminosity end. We include in 
this analysis all reverberation datasets for which (a) there is 
a reliable time lag measured for the H/3 emission line, and 



(b) there is medium V-band (F547M or F550M) HST imag- 
ing available so the host-galaxy contribution to the rest-frame 
5 100 A flux can be determined and removed. Other archival 
HST images are available for some of the objects not included 
here, but these images are not suitable for our analysis for one 
of three reasons: (1) they were taken with a different filter and 
therefore include emission lines from the galaxy and/or the 
narrow-line region, both of which would have to be corrected, 
and they would require assumptions about the unknown un- 
derlying stellar populations in the galaxy, and/or (2) the expo- 
sures are too shallow to accurately constrain the host-galaxy 
surface brightness profiles, or (3) the images are heavily sat- 
urated in the nucleus, with strong bleeding and a loss of in- 
formation at the galaxy center. We do not include reverbera- 
tion measurements of other Balmer line s in this analysis be- 
cause previous work dKaspi et al.l 120001: iBentz et al.l [2009b) 
has shown that there are differences in the mean time lags 
determined for different Balmer lines, most likely caused by 
radiative-transfer effects in the BLR clouds. We are left with 
the sample of 41 AGNs that are listed in Table Q] 

3. HST OBSERVATIONS AND DATA REDUCTION 

Nine AGNs with new H/3 time lags were imaged during 
Cycle 17 (GO- 11662, PI Bentz) with the Wide Field Cam- 
era 3 (WFC3) UVIS channel through the F547M (Stromgren 
y) filter (A e = 5447 A and AA = 650 A). This imaging setup 
allowed us to probe the continuum flux from the AGN and 
the host galaxy while avoiding strong emission lines. One or- 
bit was dedicated to each object, and each orbit was divided 
into two sets of exposures separated by a dithering maneuver 
to improve the sampling of the WFC3 point-spread function 
(PSF) and facilitate in the rejection of cosmic rays and de- 
tector artifacts (such as transient warm pixels). To maximize 
the dynamic range of the final images, each set of three expo- 
sures was graduated in time, with exposure times of approxi- 
mately 30 s, 300 s, and 690 s. We did not dither during an ex- 
posure sequence to ensure that all three images were taken at 
the same position. Most of our targets were compact enough 
to fit on a single chip of the UVIS channel, but for NGC 68 14 
we employed a larger dithering maneuver to ensure that there 
was no loss of information because of the gap between the 
chips. Details of the HST observations are given in Table [2] 

We were able to correct for saturation in the long exposures 
by making use of the linear nature of charge-coupled devices 
(CCDs). Saturated pixels in the nucleus of each galaxy were 
identified in each image by consulting the data quality frames 
from the HST pipeline. These saturated pixels were clipped 
from the image and replaced by the same pixels from a shal- 
lower, unsaturated exposure, but scaled up by the exposure- 
time ratio. Cosmic rays were cleaned with t he Laplacian 
cosm ic ray identification package L. A. Cosmic dvan Dokkuml 
2001). All of the frames for a single object were then com- 
bined with the multidrizzle task to create a distortion-free im- 
age of each AGN host galaxy. The final combined, drizzled 
images are shown in Figure Q] with the ground-based spectro- 
scopic monitoring apertures overlaid. It can easily be seen that 
the host galaxy of each AGN contributes a significant amount 
of light within the monitoring aperture. 

4. GALAXY SURFACE BRIGHTNESS DECOMPOSITIONS 

An important component of calibrating the R^lr-L rela- 
tionship is properly correcting the L measurements for the 
contribution from host-galaxy starlight. The method we em- 
ployed here is similar to that described by B entz et all (l2006al 
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2009a), where the analysis of 32 galaxies in our sample is 
reported; it relie s on using the software program Galfit 
dPeng et alj|2002l) to model the surface brightness profiles of 
the host-galaxy images. Galfit is a nonlinear least-squares 
two-dimensional image-fitting algorithm. We used the latest 
version of GALFIT (Version 3) which allows for the model- 
ing of spiral arms, rings, and irregular shapes (see iPeng et all 
2010 for a full description and various examples). 

For the surface brightness decomposition of each of the 9 
new HST host-galaxy images in this work, we emplo yed a 
tilted plane for the background sky flux and a Tiny Tim dKristl 
1993) model for the unresolved AGN. Tiny Tim models were 
generated for each specific AGN by creating a model at the 
specific detector position of each of the pointings and com- 
bining these models through multidrizzle in the same way that 
the AGN images were co mbined . Host- galaxy bulges, disks, 
and bars were all fit with Sersic (1968) profiles of the basic 
form 



S(r) = E c exp 




(1) 



where S e is the pixel surface brightness at the effective ra- 
dius r e . The Ser sic index, n, has a valu e of 1 for an exponen- 
tial disk, 4 for a ide Vaucouleursl (fl~948) profile, and 0.5 for a 
Gaussian. Bulge and bar components were modeled by allow- 
ing the Sersic index to vary with no constraints while disks 
were modeled by holding the Sersic index fixed at a value 
of 1 . Fits that resulted in bulge Sersic indices outside of the 
range ~ 0.1-6 were considered unphysical and therefore un- 
acceptable. In these cases, we required multiple PSF models 
in the center of the galaxy, offset by fractions of a pixel, to 
keep the Sersic index of the bulge from running up the maxi- 
mum value allowed by GALFIT, n = 20. A high Sersic index 
has a very peaky shape with strong wings, and can mimic a 
PSF+sky model. A runaway Sersic index in our galaxy fit- 
ting is likely because of the well-known PSF mismatch that 
can occur between Tiny Tim models and WFC3 images due to 
spacecraft "breathing" and/or jitter, but could potentially be 
caused by any marginally resolved nuclear flux from hot gas 
or star clusters. We assume here that the cause is PSF mis- 
match and ascribe all of the flux in these multiple PSF models 
(which we assume are modeling a single physical component) 
to the AGN itself, and we describe various tests of the validity 
of this assumption below. 

Galfit allows for surface brightness decompositions that 
can be as simple or complicated as the user may wish. The ul- 
timate goal of the surface brightness modeling in this project 
was to accurately remove the AGN PSF, thereby creating 
an "AGN-free" image of each host gal axy from which the 
starlight contribution could be measured. iBentz et all (12009a!) 
assume uncertainties of 0. 1 mag in the measured host-galaxy 
flux based in the range of acceptable models that could be 
found to fit an image. Here, we investigate the uncertainty 
in the best-fit models by carrying out two independent sets of 
surface brightness decompositions for each of the host-galaxy 
images in this study. 

The first set of models, which we will refer to as the "opti- 
mal" models, include multiple surface brightness components 
and make use of power-law rotation to model spiral arms, 
Fourier modes to account for "bending" of the ellipse mod- 
eling the light distribution and other asymmetric flux distri- 
butions, and truncation functions to allow for the modeling 
of rings. These models, which are shown in the bottom pan- 



els of Figures l2T[T0l are the best representations for the actual 
two-dimensional surface brightness distributions of the host 
galaxies. 

The second set of models, which we refer to as "simple" 
models, do not make use of power-law rotations, Fourier 
modes, or truncation functions. The simple models typi- 
cally require a factor of 3 fewer free parameters than the op- 
timal models and they are computationally much faster to 
run and to converge, but they less accurately represent the 
two-dimensional surface brightness profiles of the AGN host 
galaxies, as can be seen in the top panels of Figures EHTUl 
Conversely, the optimal models do a good job of reproducing 
the relative flux in each pixel in the images, but the impor- 
tance of this, other than being aesthetically pleasing, is not 
clear and the physical interpretation of each model compo- 
nent is not straightforward to determine. 

Tables l3TfTTI give the parameters determined for the "opti- 
mal" and "simple" fits for each of the 9 galaxies fit here. The 
formats of the tables are as follows: Column (1) gives the note 
for the type of fit described ("optimal" or "simple"); column 
(2) gives the component number of the fit, generally in order 
of increasing angular size and increasing angular offset from 
the center of the galaxy; and column (3) gives the description 
for the type of model component (or components in the case 
of a PSF model and tilted plane sky model). The remaining 
columns describe the various parameters of each model, with 
Column (11) listing any notes relevant to the models. We give 
a brief description of the rema ining columns bel ow, but the 
interested reader is referred to [Peng et all (2010) for further 
details of the models and their parameters employed by GAL- 
FIT. 

For the PSF models, columns (4) and (5) are the angular 
offsets in arcseconds from the center of the galaxy (defined 
as the location of the AGN PSF) in the x and y directions, 
respectively. Column (6) is the integrated magnitude of the 
PSF model. For the sky models, column (8) gives the average 
value of the sky background level in counts at the geometric 
center of the image, and columns (9) and (10) give the flux 
gradients in the x and y directions respectively. 

Sersic models are listed with columns (4) and (5) as the an- 
gular offsets in arcseconds from the center of the galaxy in 
the x and y directions, respectively. Column (6) gives the 
integrated magnitude of the Sersic component and column 
(7) lists the effective radius in arcseconds. Column (8) gives 
the Sersic index, which was held fixed at a value of 1.0 for 
exponential-disk components. Columns (9) and (10) are the 
axis ratio and the position angle of the major axis in the im- 
age. Note that images were fit at the orientation obtained dur- 
ing the observation, and the position angles listed would need 
to be corrected for the roll angle of the spacecraft to determine 
their orientation relative to north. 

For the "optimal" fits, the Sersic models were modified by 
power-law rotation, Fourier modes, and/or radial truncation 
functions. In the case of the truncation functions, where por- 
tions of the underlying model are removed, the Sersic pro- 
file is listed as "sersic3" and we report the surface bright- 
ness at the break radius (Et,) and the break radius (r\,) itself 
in columns (6) and (7), rather than the integrated magnitude 
and effective radius. 

Power-law rotations of a Sersic profile are denoted by 
"power" in column (3). Columns (5) and (6) list the inner and 
outer radii of rotation in arcseconds. Column (7) gives the ro- 
tation angle between the inner and outer radii and column (8) 
is the power-law slope, denoted as a. Columns (9) and (10) 
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are the line-of-sight inclination angle of the disk, 9i nc \ (with 
#inci = being equivalent to face-on), and the position angle 
of the rotation in the plane of the sky, 6> s k y . 

Fourier-mode modifications to Sersic profiles are denoted 
by "Fourier" in column (3). Beginning with columns (5) and 
(6), and continuing through column (10), are the modes (e.g., 
m = 1 or m = 3) and their amplitudes a m and phase angles 
(</>), where the phase angle is the relative angle between the 
Fourier mode and the position angle of the major axis of the 
parent Sersic profile. To avoid degeneracy with the axis-ratio 
parameter for the Sersic profiles, we did not make use of the 
m = 2 Fourier mode. 

Truncation functions were generally used to model rings in 
the galaxies and are denoted as "radial" in column (3). Both 
inner and outer truncations were used, with each denoted ap- 
propriately. Columns (4) and (5) give the angular offsets of 
the center of the truncation function from the center of the 
parent Sersic profile in the x and y directions, respectively. 
Column (7) gives the break radius of the truncation function, 
defined to be the radius at which the truncation function has 
a value of 99% of the flux of the untruncated Sersic model at 
that same radius. Column (8) lists the softening length, Ar so f t , 
where rb rea k ± Ar SO ft (+ for outer truncations, - for inner trun- 
cations) gives the radius at which the flux drops to 1 % of the 
untruncated Sersic model flux. Columns (9) and (10) give the 
axis ratio and position angle of the truncation function. 

Finally, the last row of each fit gives the figures of merit 
for that particular surface brightness decomposition: % 2 ; the 
number of degrees of freedom, A^of ; the number of free pa- 
rameters in the models, Nf ree ; and the reduced x 2 , xl ■ 

In addition, to test the suitability of our choice of using 
multiple Tiny Tim PSFs offset by fractions of a pixel to bet- 
ter model the AGN PSF in several objects, we carried out 
surface brightness decompositions using a variety of differ- 
ent PSF models. These included a very high S/N WFC3 im- 
age of the white dwarf EGGR 102, a bright field star from 
the image of Mrk290 (one of our targets), a fainter star in 
the field around EGGR 102, and|Mpffat ( 1969) fits to each of 
these stars. We also investigated the effect of convolving the 
image and the PSF mode l with a narrow Gaussian to ensure 
Nyquist sampling (e.g., iKim et af] 120081) . In each of these 
tests, we only allowed a single component to model the AGN 
PSF and we compare the results of the test to the results ob- 
tained using the fitting procedures described above and tab- 
ulated in Tables [3lTTlT When the image of a star was used 
as the PSF model, the difference in central host-galaxy flux 
measured from an "AGN-free" image was only ~ 1%, and it 
was only slightly higher (~ 2%) when a Moffat fit to a star 
image was used as the PSF model. Broadening the image and 
the PSF model caused the largest difference in central host- 
galaxy flux, about 7%, so while this approach has been found 
to work in the past for other HST cameras and in other situ- 
ations, it was the least successful alternative in this case. As 
expected, the fit residuals at the center of the galaxy are the 
smallest when we allow multiple Tiny Tim models to account 
for the central AGN PSF. Furthermore, Tiny Tim has the ad- 
vantage of producing PSF models with infinite S/N, therefore 
avoiding the problem of introducing additional noise into the 
"AGN-free" images from which we determine the host-galaxy 
starlight contribution. We find the same results when we also 
fit the image of the bright star in the field of Mrk 290 with 
multiple Tiny Tim PSF models — the residuals in the fit are 
decreased without adding extra noise. 



Below, we provide some notes on each of the individual 
galaxies modeled in this work. 

Mrk 142 — Mrk 142 is a late-type spiral galaxy at intermedi- 
ate inclination. It does not appear to have a bulge, but does 
seem to have a bar that manifests itself as a compact struc- 
ture with a low Sersic index (n < 1) and elongated shape. 
At the highest redshift of any of the galaxies modeled here 
(z = 0.045), the bulge may be too compact to disentangle from 
the very bright unresolved AGN. Its morphological classifica- 
tion according to these images is SBcd-SBd. The parameters 
for the optimal and simple surface brightness decompositions 
of Mrk 142 are tabulated in Table|3] and the models and resid- 
uals are displayed in Figure [2] 

SBS1116+583A — SBS 1116+583A is a relatively face-on 
barred spiral galaxy with an exponential bulge, approximately 
SBb in type. The best-fit parameters for its surface brightness 
decompositions , which are given in Table |4] include a lens 
(uniform disk, Ide Vaucouleurs|[T959l) that is more extended 
than the bulge, nearly circular, and has a very low Sersic in- 
dex (n w 0.3-0.4) in both the "optimal" and "simple" models. 
The models and residuals are displayed in Figure [3] 

Arpl51 — Arpl51 (Mrk 40) is an early-type spiral galaxy 
(SO-Sa) with a hint of remaining spiral structure and a long 
tidal tail stretching north-northwest from a recent encounter 
with a small companion galaxy at a projected angular dis- 
tance of ~ 19". The messy morphology of Arp 151 and its 
companion required multiple surface brightness components 
for an accurate fit, and we do not attempt an interpretation of 
their physical meaning here. The best-fit parameters for its 
surface brightness decompositions are listed in Table [5] and 
the models and residuals are displayed in Figure 2] 

Mrk 1310 — Mrk 1310 is a ringed spiral galaxy, approxi- 
mately Sb in type, with an apparently large number of bright 
globular clusters. There is also a faint galaxy directly south of 
Mrk 1310 that appears as an arc. The location of the galaxy 
along the line of sight to Mrk 1310 is unknown, but the dis- 
torted shape of this faint galaxy may mean that it is being 
tidally disrupted by Mrk 1310, or it may simply be a chance 
superposition.. The large angular separation and orientation 
of elongation rule out the possibility of gravitational lensing. 
Additional color information, at the minimum, will be neces- 
sary to determine where this small galaxy exists along our line 
of sight to Mrk 1310. The best- fit parameters for the surface 
brightness decompositions of Mrk 1310 are given in Table [6] 
and the models and residuals are displayed in Figure [5] 

Mrk 202 — Mrk 202 is a compact face-on spiral galaxy, ap- 
proximately Sb in type, with a bright star-forming ring. The 
best-fit parameters for its surface brightness decompositions 
are listed in Table [7] and the models and residuals are dis- 
played in Figure [6] 

NGC4253 — NGC 4253 (Mrk766) is a barred spiral galaxy of 
type SBc with a distinct nuclear spiral and a faint outer ring . It 
is also classified as a narrow-line Seyfert 1 because of its rel- 
atively narrow broad emission lines. The best-fit parameters 
for its surface brightness decompositions are given in Table[8] 
and the models and residuals are displayed in Figure [7] 

NGC 4748— NGC 4748 is a barred spiral galaxy with a nu- 
clear starbursting ring and is currently undergoing an interac- 
tion with another, slightly smaller, spiral galaxy. The best-fit 
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parameters for its surface brightness decompositions are given 
in Table|9] and the models and residuals are displayed in Fig- 
ureg] 

Mrk290 — Mrk290 is an early-type spiral galaxy (Sa-Sab) at 
a relatively low inclination to our line of sight. The bright 
source to the southeast appears to be a star in our own galaxy. 
The best-fit parameters for its surface brightness decomposi- 
tions are listed in Table [10] and the models and residuals are 
displayed in Figure [9] 

NGC6814— NGC6814 is a beautiful, face-on, moderately 
barred spiral galaxy at a fairly low redshift of 0.0052 (Dl « 
20 Mpc), making it one of the nearest broad-lined AGNs in the 
local universe and in our sample. The best-fit parameters for 
its surface brightness decompositions are given in Table QT| 
and the models and residuals are displayed in Figure [TOl 

5. AGN FLUXES 
5.1. Starlight Measurements and AGN Flux Recovery 

The host-galaxy starlight contribution to the 5 100 A spec- 
troscopic flux was determined by first measuring the yield 
of electrons within a rectangular aperture, with dimensions 
and orientation matching that of the ground-based moni- 
toring campaign, centered on the nucleus of the galaxy 
in the PSF- and sky-subtracted Hubble Space Telescope 
(HST) image. The exposure time and inverse sensitivity 
for each image {HST keyword photflam, having units of 
ergs cm" 2 A" 1 electron" 1 ) were utilized to recover the incident 
photon flux from the yield of electrons. All photflam values 
were taken from the most recent recalibration of the appropri- 
ate dataset through the HST pipeline as of 2012 June 12. For 
the ACS images, the photflam values are somewh at different 
from those previously used by iBentz et ail (12009a!) because 
they have been updated to account for th e loss of sensitivity of 
the High Resolution Channel over time (Bo hlin et al.112.01 II) . 

Once the host-galaxy flux through the HST system had been 
determined, a small color correction was necessary to account 
for the difference between rest-frame 5 100 A and the pivot 
wavelength^ of the filter. To determine the c olor correction, 
a bulge template spectrum ( Ki nney et al.l 1996) was redshifted 
and reddened by the appropriate amounts to approximate the 
central host galaxy of each AGN. Only Galactic extinction 
was included in the reddening, and the values used are slightly 
differen t from previous values emp loyed by |Bentz et al.l 
(2009a) because they are based on the Schlafly & Finkbeiner 
(120111) recalibration of the ISchlegel et al.l (119981) dust map. 
The Galactic extinction values are smaller by a few hun- 
dredths of a magnitude for all of our sample (median differ- 
ence of -0.024 mag) except for 3C 120 where the new extinc- 
tion value is 0.2 mag smaller than before. The ratio of the 
5100 A flux to the flux through the HST filter (which we asso- 
ciate with the pivot wavelength) was estimated from the red- 
shifted and reddened spectrum using synphot, and is listed 
in Table [T2l The final derived host-galaxy flux contributions 
to the ground-based spectroscopic continuum flux are given 
in Table [T2l These values were subtracted from the absolute 
calibrations of the mean continuum fluxes during the monitor- 
ing campaigns to recover the mean AGN fluxes at rest-frame 
5 100 A, from which the AGN luminosities were determined. 

18 A measure of the effective wavelength of a filter that is independent of 
the source spectral energy distribution (Tokunaea & Vacca 2005). 



We discuss the effects of ground-based seeing and modeling 
uncertainties, among others, below. 

5.2. Uncertainties 

The uncertainty in the recovered AGN flux is a combina- 
tion of the mean measurement uncertainty in the continuum 
flux from the reverberation campaign and the uncertainty in 
the host-galaxy contribution to the continuum flux. The for- 
mer is a small component, ranging from 1-5% for the mea- 
surements included here. This is due to the requirement for 
reverberation-mapping campaigns to achieve a high S/N per 
pixel in the continuum flux of each individual spectrum ac- 
quired throughout the campaign (typically S/N « 30-100) in 
order to measure the few-percent variations that evidence the 
reverberation signal. The latter contribution to the AGN flux 
uncertainty was determined by adding in quadrature the un- 
certainty in the starlight flux from the modeling and the un- 
certainty in the starlight flux from ground-based seeing effects 
that would be in place during a reverberation-mapping cam- 
paign. 

Modeling Uncertainties — The uncertainty from surface 
brightness modeling was determined by comparing the 
starlight measurement derived for each object from the sim- 
ple fit and the best fit detailed in the previous section. The 
addition of Fourier modes and power-law rotation, in gen- 
eral, changed the starlight measurement by 0.04% (Mrk202) 
to 8% (Mrk766), with a median difference of 3%. Because 
such comparisons are time-intensive and computationally de- 
manding, we have not carried them out for all 41 galaxies in 
the sample. Instead, we adopt a conservative estimate of 5% 
uncertainty for the host-galaxy contribution for all compact 
galaxies, where the field of view of the HST image contains 
a large fraction of pixels that consist of empty sky (e.g., the 
Markarian objects and the PG quasars). For extended galaxies 
that fill the field of view of the HST camera with which they 
were observed (e.g., the NGC galaxies with the ACS HRC), 
we adopt a 10% uncertainty in the host-galaxy contribution 
due to the greater uncertainty in the determination of the back- 
ground sky level during the modeling process. 

Seeing Effects — Optical reverberation-mapping campaigns 
are generally carried out from the ground and thus have to 
contend with variable seeing from night to night throughout a 
campaign. The effects of slit losses and variable seeing on the 
measured AGN flux are minimized by using a wide spectro- 
scopic slit (4"-5") and by carrying out an internal calibration 
of all the spectra obtained for an object by utilizing the non- 
variable [O III] AA 4959, 5007 doublet. Nevertheless, seeing 
redistributes the galaxy flux as well and can cause the starlight 
measurements from diffraction-limited HST images to dif- 
fer from the contribution obtained through the ground-based 
setup under typical seeing conditions. Reverberation cam- 
paigns generally scale the final spectra to the [O III] flux mea- 
sured on photometric nights, with a typical seeing of ~ 1". 
To investigate the effect of seeing on the derived host-galaxy 
flux, we took the "AGN-free" images of NGC 5548 (an ex- 
tended galaxy) and SBS 1116+583A (a compact galaxy) and 
created a simulated ground-based image of each by smear- 
ing with a 1" FWHM Gaussian. The starlight measurements 
were then made in the same way from the simulated ground- 
based images as they were from the diffraction-limited im- 
ages. The difference in measured starlight flux was negligible 
for NGC 5548, only 2%, but was 8% for SBS 1 1 16+583A. 
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Based on these results, we adopt an average 5% uncertainty 
for the host-galaxy contribution for each object in the sample 
due to ground-based seeing effects. 

Background Determination — Finally, we have also considered 
the effect of background subtraction during the spectral re- 
ductions on the host-galaxy flux. For the extended galaxy 
NGC 5548, we measured the host-galaxy flux in the "back- 
ground" regions on either side of the extraction region. The 
average of these background regions was treated as "sky" flux 
and subtracted from the flux within the extraction region. The 
difference in host-galaxy flux was found to be only 2% even 
though NGC 5548 is a bright extended galaxy. For the more 
compact galaxies in our sample, the effect would be even less. 
Therefore, we consider the effect of background-subtraction 
regions during spectral reductions to be negligible. 

6. DISTANCES AND AGN LUMINOSITIES 

By far, the largest contribution to the uncertainty of the 
AGN luminosities is from the uncertain distance to each 
AGN. Only five of the 41 AGNs in this study have distance 
measurements independent of their redshifts — NGC 3227, 
NGC 3783, NGC 4051, NGC4151, and NGC 4593 — while 
for the remaining 36 we estimate the distance from the red- 
shift of the AGN. The distance measurements for the five 
aforementioned objects generally come from an average of 
the distance moduli for galaxies within the same group and 
were gen erated as part of a study of the "local" velocity 
anomaly dTullv e t al. 2008)); the y were retrieved f rom the Ex- 
tragalactic Distance Database dTullv et al.l 12009b . They are 
calibrated to the same zeropoint as the HST Key Project 
dFreedman et al.ll2001h . which found Ho = 72kms _I Mpc -1 , 
and which we have adopted throughout this work. 

In general, we find that the uncertainties in the distances are 
underestimated for these individual sources in the Extragalac- 
tic Distance Database. For each of these five objects, we next 
give a brief discussion of the available distance measurements 
and their apparent quality, as well as the distances we adopt. 

NGC 3227 — There are seven galaxies in the same group as 
NGC 3227, with distance measurements ranging from 18- 
34 Mpc. Fortuitously, however, NGC 3227 is currently inter- 
acting with the early-type galaxy NGC 3226. The distance 
to NGC 3226 from the sur face brightness flu ctuations (SBF) 
method is 23.5 ± 2.4Mpc dTonrv et al .1120011) , which is 10% 
less than the group-averaged distance estimate to NGC 3227 
of 26.4 ± 1.6 Mpc. We adopt the distance measurement of 
NGC 3226 as the distance to NGC 3227. 

NGC3783 — The three galaxies in the group to which 
NGC 3783 belongs have measured distances ranging from 
20 to 28 Mpc, leading to a group-averaged distance esti- 
mate of 25.1 ±2.9 Mpc for NGC 3783. Based on its red- 
shift of 0.00973, however, NGC 3783 is estimated to lie at 
a distance of 41 Mpc. This is a difference of nearly 50% 
in distance that translates into a factor of almost 3 differ- 
ence in predicted luminosity. With a recessional velocity of 
29 17 km s" 1 , NGC 3783 would generally be expected to have 
peculiar velocities affecting its perceived recessi onal veloc- 
ity at only the ~ 10% leve l (~ 300 km s" 1 , e.g., Masters et al. 
120061 iBahcall & Oh|[T996l) . a severe underestimate given the 
50% discrepancy between the group-averaged distance esti- 
mate and the estimate based on redshift. We adopt the group- 
averaged distance of 25.1 Mpc, with an uncertainty of 20% 
(5.0 Mpc), for NGC 3783. 



NGC 4051— NGC 4051 is one of 64 galaxies identified as 
belonging to the same group. The Tully-Fisher distance to 
NGC 4051 is quoted as 12.2 Mpc, but it does not appear to 
have been corrected for the contribution of the AGN to the 
total galaxy luminosity. The AGN contribution would appear 
to make the galaxy brighter, and it would therefore seem to 
be nearer than it actually is. The group-averaged distance 
of 17.1 ±0.8 Mpc includes individual galaxy distances rang- 
ing from 10-30 Mpc. Within the group of 64 galaxies, there 
is one galaxy with a Cepheid distance and eight early-type 
galaxies with distances from SBFs. The distances for these 
nine galaxies, which are expected to be more accurate on an 
individual basis than the distances to the other 55 galaxies in 
the group, span a smaller range of 10-21 Mpc, with a median 
value of 14.3 Mpc that is fairly consistent with the average 
distance found for the full group of 64 galaxies. Based on 
our limited information regarding the location of NGC 405 1 
within its group, we adopt the group-averaged distance of 

17.1 Mpc for NGC 405 1 based on all the galaxies in the same 
group, with an uncertainty of 3.4 Mpc (20%). 

NGC 41 51 — There are only four galaxies contributing to the 
group-averaged distance for NGC 4151, and their individual 
distances are estimated to range from 3.9 Mpc to 34.0 Mpc 
based on the iTullv & Fisheri (11977b line width-luminosity 
correlation, with a final distance for NGC 4151 quoted as 

1 1 .2 ± 1 . 1 Mpc. The object with the smallest estimated dis- 
tance of 3.9 Mpc is NGC 4151 itself, but the total galaxy lu- 
minosity does not appear to have been corrected for the enor- 
mous contribution from the AGN. It therefore appears that 
the distance of 3.9 Mpc is a gross underestimate caused by 
neglecting the AGN contribution to the total galaxy lumi- 
nosity, causing the galaxy to appear brighter (and therefore 
nearer) than it actually is. We have recalculated the group- 
averaged distance while excluding the likely erroneous dis- 
tance of 3.9 Mpc and adopt a distance of 16.6 ±3.3 Mpc for 
NGC 4151. 

NGC 4593 — Only two galaxies contribute to the group- 
averaged distance for NGC 4593. They have individual dis- 
tance measurements of 33 Mpc and 43 Mpc. We adopt the 
distance estimated by averaging their distance moduli and an 
uncertainty of 20% (37.3 ±7.5 Mpc) for NGC 4593. This is 
consistent with the distance of ~ 39 Mpc expected from the 
redshift of NGC 4593 and assuming that NGC 4593 has zero 
peculiar velocity. The 20% uncertainty in the distance that we 
have assumed for NGC 4593 (and, indeed, several of the other 
distance estimates above) may be an underestimate of the true 
discrepancy between the actual distance to the source and our 
estimate of the distance. 

For the other 36 AGNs in the sample, we have no choice at 
this time but to estimate their distances from their measured 
redshifts. Because of this, peculiar velocities can introduce 
a large uncertainty into these distance estimates. To further 
complicate the issue, peculiar velocities are highly direction- 
dependent (e.g., the "Finger of God" effect) and will be ran- 
domly oriented relative to our line of sight for the AGN host 
galaxies in this sample. Our lack of additional distance infor- 
mation for the vast majority of the AGNs in the reverberation 
sample leads us to conservatively estimate that peculiar veloc- 
ities affect the galaxy recession velocities at an average level 
of - 500 kms" 1 , or ~ 17% for z = 0.01. We caution that this 
may still be a significant underestimate of the accuracy of our 
assumed distances for some individual galaxies, as it would 
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be in the above case of NGC 3783 if we had no information 
beyond the galaxy's redshift. 

Clearly, there is a desperate need for accurate distance mea- 
surements to the AGN host galaxies in the reverberation- 
mapping sample. The Tully-Fisher method has been shown 
to be accurate to ~ 20% for individual galaxies and can reach 
spiral galaxies out to z ~ 0.1, but will require extra care for 
these galaxies to ensure removal of the AGN contribution to 
the total galaxy luminosity. Furthermore, there are a handful 
of galaxies in the sample that are within reach of the ~ 30 Mpc 
limit for Cepheid observations with HST. We have an ap- 
proved Cycle 20 program to obtain a Cepheid-based distance 
measurement to the face-on spiral galaxy NGC 68 14 (GO- 
12961, PI Bentz). NGC4151, in particular, is another galaxy 
with a very large distance uncertainty that would benefit from, 
and be within the reach of, an HST Cepheid program. 

Finally, we note that we do not attempt to correct for in- 
ternal reddening from the AGN host ga laxy. Previous studies 
dBentz et alj2009allDennev"eT al. 2010) have shown that such 
corrections, for the few objects where they are possible, are 
fairly small relative to the large distance uncertainties we have 
described above. In the case of the reddened AGN NGC 3227, 
for example, the reddening curve derived by Crensha w et al.l 
(120011) gives an extinction of 0.26 dex in luminosity at 5 100 A. 

Table [T3l lists the 5100 A luminosities we have determined 
for each of the datasets in our sample using the distances dis- 
cussed above. We also give the corresponding broad H/3 time 
delays, which we take to be the average radius of the H/3- 
emitting BLR. 

7. THE RADIUS-LUMINOSITY RELATIONSHIP 

Based on previous work dKaspi et all HOOOI 120051: 
IBentz et alj|2006all2009ah . we expect the form of the Rblr-L 
relationship to be a power law. We parametrize the Rblr-L 
relationship here as 

log(flBLR/l It-day) = tf+alog(AL A /10 44 ergss" 1 ). (2) 

To determine the best fit to the Rmr -L relations hip, we 
employed the LINMIX_ERR algorithm (iKellyl l2007h . which 
takes a Bayesian approach to linear regression with mea- 
surement errors in both c oordinates a nd a component of in- 
trinsic, random scatter. Kelly] (12007b carried out extensive 
tests of the consistency between LINMIX_ERR an d the com- 
monly used algorithm s FITEXY dPress etal.ll 992!) and BCES 
(Akritas & Bers hadyiri996l) . finding that the best fits deter- 
mined by all of the algorithms were generally consistent, but 
that even in cases of large scatter or poorly constrained mea- 
surements, LINMIX_ERR always derived a fit that was consis- 
tent with the known parent population from which the mea- 
surements were sampled. We fit the Rblr-L relationship with 
FITEXY and BCES and found that all the algorithms pro- 
vided consistent results, as we expected. We report the best- 
fit parameters determined by the LINMIX_ERR algorithm in 
Table [p 

The issue of dealing with multiple measurements for a sin- 
gle object when fitting the Rblr-L relationship is not straight- 
forward. On the one hand, if an individual AGN moves along 
its own Rblr-L relationship that is parallel to the Rblr-L rela- 
tionship investigated here, then there is no real difference be- 
tween multiple measurements of a single object versus mea- 
surements of many different objects. In this case, all measure- 
ments should be given equal weight in the regression analysis, 
regardless of which AGN they "belong" to. If, on the other 



hand, individual AGNs have individual Rblr-L relationships 
that are oriented at some other slope relative to the Rblr-L 
relationship for the population, then each AGN should only 
be allowed to contribute a single measurement to the regres- 
sion analysis. This issue has been exa mined for th e AGN 
NGC 5548 by iPeterson et all dlool and IBentz et all J2007T) . 
and both studies found that the slope of the optical Rblr-L 
relationship for NGC 5548 alone is steeper than the global re- 
lationship. Preliminary work by Kilerci Eser et al. (in prepa- 
ration) is also finding that the slope of the optical Rblr-L 
relationship is steeper for an individual object with multiple 
measurements and that the scatter introduced into the global 
relationship appears to be modest, but these results are neces- 
sarily based on only the small number of targets with multiple 
time lag measurements. We have also, therefore, investigated 
the effect of randomly choosing only a single measurement 
to represent each object in the sample, and the best-fit Rblr- 
L relationship is consistent within the uncertainties with the 
best-fit relationship derived with all the measurements for ev- 
ery AGN included. 

Figure QT| displays the current version of the H/3 Rblr-L 
relationship. In the top-left panel, all of the individual mea- 
surements are plotted. The new measurements that have been 
added to the Rblr-L relationship in this work are shown as 
open circles and preferentially populate the low-luminosity 
end of the relationship. 

7.1. Notes on Individual Objects 

Mrk 142 — The narrow-line Seyfert 1 galaxy Mrk 142 appears 
to be a significant outlier in the Rblr-L relationship. With an 
AGN luminosity of 3.5 x 10 43 ergs s" 1 , it has a predicted H/3 
time lag of ~ 20 days, compared to its observed time lag of 
~ 3 days. In spection of t he con tinuum and H/3 light curves 
presented by IBentz et ail (12009b!) shows a lack of strong fea- 
tures, such that no apparent lag can be detected by eye. Fur- 
thermore, the cross-correlation function for Mrk 142 shows 
the lowest significance for the objects in the LAMP sample 
with reported lag detections. We suggest that the reverbera- 
tion experiment for Mrk 142 should be repeated in an effort 
to detect a more significant time lag. Repeating the rever- 
beration experiment will also allow for the confirmation or 
contradiction of the outlier status of Mrk 142. While it is 
possible that Mrk 142 is truly an outlier, it is worth noting 
that other narrow-line Seyfert Is in the reverberation sam- 
ple (NGC 4051, NGC 4253, NGC 4748) lie extremely close 
to their expected locations, and well within the sample scat- 
ter. 

Arpl51 — The tidally distorted galaxy Arpl51 was one of 
the most variable objects in the 2008 LAMP campaign, with 
a well-determined time lag of 4.0 ± 0.5 days. However, there 
appears to be a problem with the flux calibration for Arp 151. 
The flux calibration for the LAMP sample of AGNs was de- 
termined by comparing the observed [O III] A5007 A flux for 
NGC 5548 to the known [O III] flux from many years of spec- 
trophotometric monitoring. Because of the unstable nature of 
weather and the need for many observations over a long pe- 
riod of time, reverberation datasets rely on the narrow [O III] 
emission lines as internal calibration sources for the many 
nonphotometric nights on which data is obtained. The emis- 
sion lines do not vary on the timescales probed in a single 
monitoring campaign, and they provide the final multiplica- 
tive flux calibration factor in the spectral pipeline. 
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The [O III] flux of NGC5548 is well-known because of 
the many years that this AGN has been monitored and it was 
the only well-studied AGN included in the 2008 LAMP cam- 
paign. Looking at the measured [O III] fluxes for every night 
during the 2008 campaign, there seemed to be only a single 
night among the 64 nights of the campaign where the weather 
at Lick Observatory was steadily photometric, providing an 
accurate [O III] flux measurement for NGC5548. Unfortu- 
nately, it appears that this may not have actually been the 
case, at least during the observations of Arp 151. The problem 
arises from the fact that any reasonable flux for the AGN in the 
ffi>r image of Arp 151 is more than the total continuum flux 
of AGN+galaxy derived from the ground-based spectroscopy. 

Does this mean that the process of estimating host-galaxy 
flux from HST imaging is flawed? Probably not. The clue 
comes from reviewing the final spectroscopic flux calibrations 
for the LAMP AGNs. The flux in the mean spectrum for each 
object required an increase by a multiplicative factor to match 
the measured [O III] flux from the single photometric night. 
In the case of Arp 151 (and only Ar,151), however, and only 
Arp 151, the [O III] flux measured from the supposed pho- 
tometric night of observations is less than the mean [O III] 
flux from all of the observations. Reducing the mean flux by 
the derived multiplicative factor results in a mean continuum 
flux that is too small and thus produces a negative AGN flux 
when the host-galaxy starlight correction is applied. If we 
instead set the multiplicative factor equal to 1, the problem 
disappears. Based on the multiplicative correction factors for 
all of the other objects in the LAMP sample, we should expect 
that this factor is > 1 for Arp 151 . 

Indeed, in reviewing the weather logs from the 2008 LAMP 
campaign and focusing only on nights logged as possibly 
photometric and with [O III] fluxes that are consistent with 
each other, we derive a flux of 0.76 x 10" 13 ergs s" 1 cm" 2 . 
This is 56% larger than the originally derived [O III] flux of 
0.49 x 10~ 13 ergs s" 1 cm" 2 . Furthermore, it is relatively con- 
sistent with the [O III] flux of 0.83 x 10" 13 ergs s" 1 cm" 2 de- 
rived fr om potentially ph otometric nights during the 2011 
LAMP (iBarth et al.ll201 ll) monitoring of Arp 151, given the 
slightly larger extraction width adopted during that campaign 
(10"3 compared to the 9 '.'A extraction width used in the 2008 
campaign) and the overall uncertainty in the flux calibration 
for this single object. 

We therefore adopt a flux correction factor of 1.56 as de- 
rived above, implying that the mean continuum flux den- 
sity at rest-frame 5100A for Arp 151 was (1.835 ±0.079) x 
10" 15 ergs s" 1 cm" 2 A during the LAMP 2008 campaign. We 
include this corrected measurement in the values tabulated in 
Table E] 

PG 2130+099 — There have been several reverberation exper- 
iments targeting PG 2 130+099, yet there remain s some am- 
biguity as to the accuracy of the reported lags. Kas pi et al.1 
(2000) determined a tim e lag of ~ 168 days, but a reanaly- 
sis bv lGrier et al.l (1200 8 ) found evidence for aliasing based on 
seasonal gaps in the light curve. Analysis of individual sea- 
sonal light curves gave much shorter time lags (16^-4 days). 
Two new reverbe ration experiment s have led to re ported time 
lags o f 23 days dGrier et al.l 120081) and 10 days (iGrier et al.l 
l2012al) . The latter experiment is the most recent and has the 
best time sampling for this object to date: Af me d = 0.5 days 
for the continuum and Af me< i = 1.0 days for the emission line. 
With a monitoring baseline of ~ 120 days, however, and a 



predicted time lag of ~ 40 days based on the luminosity of 
PG 2 130+099, in retrospect the experiment is uncomfortably 
close to the minimum time baseline reco mmended for rever- 
beration experiments dHorne e t al. 2004). The most obvious 
feature in the continuum light curv e occurs at a helio centric 
Julian day of HJD-2450000 » 55 10 (IGrier et al.l2012al Figure 
2) and would be expected to be echoed in the H/3 light curve 
at HJD rs 5550, right where th e campaign abruptly ends. Fur- 
thermore, IGrier et al.l (120 12bl) were able to reconstruct a map 
of the time-delay response as a function of velocity across 
the emission-line p rofile using the light curves presented by 
IGrier et ail (l2012al) . From this analysis, they determine that 
the H/3 time lag associated with this monitoring dataset is 
likely to be ~ 31 days, not 13 days. It is therefore debatable 
whether PG 2 130+099 is truly an outlier, and so we recom- 
mend that yet another reverberation experiment be dedicated 
to this object. The tightness of the R^lr-L relationship to 
date implies that the discovery of a true outlier may well give 
important clues about detailed AGN physics deep within the 
potential well of the central black hole. 

Because of the plausibly erroneous nature of the BLR 
radius measurements for Mrk 142 and PG 2 130+099, we 
have also carried out fits to the Rblr-L relationship with 
Mrk 142 excluded and with an adopted l ag of 3 1 ± 4 days fo r 
PG2130+099 (based on the analysis of IGrier et all 1201 2bl) . 
With these two changes, the slope of the relationship is 
slightly increased, but still consistent within the uncertain- 
ties (see Table fl4l. We have also considered the fit with both 
Mrk 142 and PG 2130+099 excluded (see TableE} and the fit 
is still consistent within the uncertainties. 

Finally, we have also investigated the effect of correcting 
for internal extinction in the one AGN where we observe a 
large reddening and we have available an appropriate red- 
de ning curve, NGC 3227 . As previously mentioned, work 
by lCrenshaw et al.l ((2001) leads to an extinction correction of 
0.26 dex at 5100 A for NGC 3227. We have applied this cor- 
rection to the luminosity measurement for NGC 3227, both as 
a member of the full sample of 7 1 measurements, and with the 
different treatments of Mrk 142 and PG 2 130+099 described 
above. The results are given in Table fl4l and again, in both 
of these cases, the changes to the best-fit solution are mini- 
mal. NGC 3227 is known to be one of the most heavily red- 
dened objects in our sample, thus, the internal extinction cor- 
rection will be much smaller for the rest of the AGNs. The 
top right panel of Figure QT| displays the Rq^r-L relation- 
ship with Mrk 142 excluded, an adopted lag of 31 ±4 days 
for PG 2130+099, and a reddening correction for NGC 3227. 



8. DISCUSSION 

The best-fit slope of a = 0.53 3^n n 3 ^ is consistent with th e 
analyses previously presented by iBentz et al.l (l2006al l2009al) . 
There appears to be no difference in the relationship at the 
high-luminosity and low-luminosity ends, with no evidence 
for a turnover at low luminosities. In fact, the Rblr-L rela- 
tionship appears to be remarkably consistent over four orders 
of magnitude in luminosity among these AGNs. Furthermore, 
the relationship is remarkably consistent with the expectation 
from simple photoi onization argum ents. Specifically, it was 
first pointed out by Davi dson] (119721) that we can define the 
ionization parameter of a BLR cloud as 



U = 



Q(H) 
ATrR 2 cn e 



(3) 
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where R is the distance from the central source, c is the speed 
of light, n e is the electron number density, and 



hv 



-dv 



(4) 



is the flux of hydrogen ionizing photons emitted by the central 
source. Under the assumptions that the ionization parameters 
and particle densities are about the same for all AGNs, one 
finds that 

R oc <2(H) 1/2 , (5) 

so that the radius at which a particular emission line is most 
likely to be emitted is a simple function of the intensity of the 
ionizing flux. Further assuming that the ionizing continuum 
shape is not a function of luminosity, such that L oc Q(H), we 
expect 

R(xL 1/2 . (6) 

The above arguments certainly gloss over many of the finer 
details of BLR photoionization physics, but this seems to mat- 
ter little as the observed relationship matches this simplistic 
expectation quite well. This particular prediction of photoion- 
ization physics, namely that the size of the BLR should scale 
with the luminosity of the central source, was sought in the 
very early days of reverberation-mappin g experiments, when 
the fir st BLR sizes were being measured (Koratkar & Gaskell 
119911) . It took another decade, however, for the BLR mea- 
surements to span a sufficiently large dynamic range so that 
the rela tionship was clear ly detected, despite the large initial 
scatter dKaspi et al.ll2000t) . 

More recent work on the physical basis for an Rblr-L re- 
lationship, spurred on by the initial and continuing successes 
of the reverberation-mapping method includi ng near-infrared 
(IR) photometric reverberation mapping (Suganuma et al. 
2006), has focused on the role of dust and the dust sublima- 
tion radius in setting the size of the BLR. Th e importance of 
dust was first noted by Netz er & Laorl (119931) and has been 
analyzed more recently bv lGoad et al.l (120121) . among others. 
The upshot of many of these models is that the outer edge of 
the BLR is bounded by the dust sublimation radius, perhaps 
coincident with the in ner edge of the du sty torus-like structure 
of the unified model d Antonuccill 19931) . Outside the dust sub- 
limation radius, the line emission from the dusty gas is sup- 
pressed by a large factor because the dust grains absorb many 
of the incoming ionizing photons as well as the emitted line 
photons, effectively creating an outer edge for the BLR. A nat- 
ural consequence of a central ionizing source with a variable 
flux is that the dust sublimation radius will respond to these 
flux variations. An increase in ionizing photons will destroy 
many dust grains and increase the dust sublimation radius, 
whereas a decrease in ionizing flux will allow more grains 
to condense or migrate in and decrease the dust sublimation 
radius. While the basic physical motivation for a Rblr-L re- 
lationship seems to be understood, there a re many details that 
are currently unknown. iGoad et ail (120121) provide a compre- 
hensive overview of the state of photoionization models and 
the agreement (or lack, thereof) with observations. 

The form of the relationship appears to be fairly well- 
determined at this point and has not changed significantly 
with the updates and additions included here. The regression 
results are also consistent with the results of a microlensing 
analysis of the BLR in lensed quasars, where the magnifica- 
tion amplitude is dependent on the size of the emission region, 



an independent method that is not subject to the same un- 
certai nties involved in reverberation mapping dGuerras et al.l 
120121) . Furthermore, the scatter about the Rblr-L relationship 
is now quite low. The LINMIX_ERR routine provides an es- 
timate of the scatter about the relationship of 0. 19 ± 0.02 dex 
(about 56%). We plot the residuals of the measured BLR radii 
to the estimated BLR radii derived from the results of the LIN- 
MIX_ERR routine fit in the bottom left panel of Figure QT| 
The residuals are approximately normally distributed, as can 
be seen by comparison to an overplotted Gaussian function 
with a = 0.19 dex (dotted line). The scatter about the rela- 
tionship is often called the "intrinsic" scatter, but in this case 
it is actually a combination of the real intrinsic scatter and 
variance from inaccurate or biased measurements. Thus, the 
intrinsic scatter in the Rblr-L relationship is likely to be less 
than the 0.19 ±0.02 dex we find here. Indeed, if we omit the 
two most suspicious measurements in the sample, those of 
Mrkl42 and PG 2130+099 as discussed above, we find that 
the scatter drops to 0. 13 dex. The residuals with the exclusion 
of Mrk 142, the adopted lag of 3 1 ± 4 days for PG 2 1 30+099, 
and the reddening correction for NGC 3227 are plotted in the 
bottom-right panel of Figure QT| This is significantly lower 
than the typical scatter in black hole scaling relationships, 
such as the 0.4 dex scatter (or larger) of the M BH - a* relation- 
ship (iGultekin et al.ll2009b iPark et al.ll2012l: iMcConnell & Mal 
12011) . 

Such low scatter in the Rblr-L relationship and the poten- 
tial to decrease it even further with accurate distances and ad- 
ditional reverberatio n campaigns seems t o lend support to the 
recent arguments of Wat son et al.l (1201 ll) that the Rblr-L re- 
lationship could potentially be used to turn any AGN with a 
well-determined H/3 time lag into a standardizable candle for 
use in cosmological studies. While there is still work to be 
done before the Rblr-L relationship can match the 0.05 dex 
scatte r in the Type la supernova Hubble diagram res iduals 
(e.g.. ISilverman et al.ll20ll iGaneshalingam et al.ll2013h . the 
main weakness at the moment is the lack of accurate dis- 
tance measurements to tie the current sample of reverberation- 
mapped AGNs onto the well-established nearby distance lad- 
der. However, once this has been rectified, the radius and flux 
of any AGN could be measured and compared to the expected 
luminosity distance derived from the Rblr-L relationship. In 
principle, this is possible out to z ~ 4 for near-IR spectroscopy 
of the H/3 emission line, far beyond the current reach of Type 
la supernovae. The high luminosities of such quasars, how- 
ever, mean that the observed H/3 time lags would be on the 
order of decades. Because the BLR is ionization stratified, 
emissions lines with a higher ionization potential than H/3, 
such as C IV or He II, have rest-frame time delays that are 
factors of a few smaller than those of H/3. If separate Rblr-L 
relationships for these lines could be define d and calibrated as 
accurately as we have achieved for H/3 (see Kas pi et ai]l2007l 
for current progress on the C IV Rblr-L relationship), then ex- 
ploring cosmology with reverberation experiments would be- 
come much more feasible for objects with z > 2, where the ex- 
pected observed time delays suffer heavily from time-dilation 
effects. The ability to probe out to such high redshifts would 
provide observational constraints on the evolution of the dark 
energy equation-of-state parameter, as well as on alternative 
theories of gravity (King et al. in preparation). 

9. FUTURE WORK 

The most pressing deficiency in the H/3 Rblr-L relation- 
ship is the current lack of accurate distance measurements 
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to the AGN host galaxies. The uncertainty in the distances 
at present provides the single largest source of un certainty in 
the AG N luminosity measurements, especially as Tul lv et all 
(2008) have shown that peculiar velocities may still be impor- 
tant (> 10%) even beyond 50Mpc. We are working to obtain 
distances based on the Tully-Fisher method for the intermedi- 
ate inclination spirals in our sample with new HI spectroscopy 
and near-infrared imaging. We will also obtain a Cepheid- 
based distance for NGC6814 (HST GO- 12961, PI Bentz), and 
will explore additional distance indicators (such as the globu- 
lar cluster luminosity function and the planetary nebula lumi- 
nosity function) for the remaining AGN host galaxies in the 
sample. New reverberation campaigns should be dedicated to 
studying the two most likely outliers in the R^lr-L relation- 
ship — PG 2 130+099 and Mrk 142 — to determine whether 
they are truly outliers. Additionally, the difficulties with the 
flux calibration of Arp 151 suggest that it would be a useful 
exercise to acquire new spectra of all the objects in the sample, 
taken under stable photometric conditions and with a uniform 
setup. 

There is also much work to be done to determine indepen- 
dent Rblr—L relationships for C IV and other emission lines 
for use in estimating black hole masses at higher redshifts, 
where the H/3 emission line has redshifted out of the observed 
optical bandpass. All of the derived Rrlr-L and black hole 
mass scaling relationships for all other emission lines in the 
literature currently rely on the H/3 R^lr-L relationship. It is 
therefore critical that we build up independent R^lr-L rela- 
tionships for these other commonly utilized emission lines. 
While Mg II and C IV are frequently employed for black hole 
mass estimates at z > 0.5, there are only a handful of exist- 
ing C IV reverberation results (see Kaspi et al. 2007), with 
the vast majority centered around a very small range in lu- 
minosities, and in the case of Mg II, only a single object has a 
measured reverberation time lag. Furthermore, the C IV time 
lags are generally deduced from International Ultraviolet Ex- 
plorer (IUE) spectra that were obtained every few days, and 
so have relatively poor tempo ral sampling comp ared to what 
is typically achieved for H/3 dPeterson et al.l20 04). The higher 
ionization state of C IV means that we expect its time lag to be 
a factor of 2-3 shorter than that of H/3. Future C IV reverbera- 
tion experiments of low- to moderate-luminosity AGNs in the 
nearby universe will require daily sampling or better in order 
to measure C IV time delays to the same level of significance 
that is now typically achieved for H/3. 

Finally, although this empirical relationship is well mea- 
sured and there seems to be some theoretical understanding 
behind it, the details of the photoionization physics, as well 
as the geometry and kinematics of the gas, are not well under- 
stood at this time. There is significant room for improvement 
in our physical understanding of AGN BLRs. 

10. SUMMARY 

We have carried out an imaging program with HST to pro- 
vide starlight corrections to the luminosities of 9 AGNs with 
H/3 radius measurements. We have fully updated and revised 
the calibration sample for the Rblr-L relationship, including 
20 new H/3 BLR measurements from recent reverberation- 
mapping campaigns, and we have reinvestigated the form of 
the relationship. We find a best fit of log(7?BLR/Ht-day) = 
1.527^1 +0.533!^log(AWl0 44 L Q ). This is consistent 
with a slope of 0.5 and with previous work that included 
starlight corrections to the AGN luminosity measurements. 
After including the additions and updates, the single largest 



source of uncertainty comes from the highly uncertain dis- 
tances to the AGNs in the sample. The low scatter in the rela- 
tionship (0.19 ± 0.02 dex) and the potential to further reduce 
the scatter, with no clear outliers, support the proposed use of 
the Rslr-L relationship to probe the matter and energy con- 
tent of the Universe out to z ~ 0.6 with optical measurements 
of the H/3 emission line. Pushing the H/3 observations into 
the near-IR would allow the relationship to probe quasars out 
to z ~ 4, beyond the reach of Type la supemovae and into a 
new interesting regime for tests of the predictions of different 
cosmological models. 
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FIG . 1 . — HST WFC3 F547M images of the AGN host galaxies, displayed with an inverted logarithmic stretch. The black rectangles show the geometry and 
orientation of each ground-based spectroscopic monitoring aperture. The size of the region displayed is l'xl', except for NGC6814 which is displayed in a 
2' x2' box. For all images, north is up and east is to the left. 



FIG. 2. — Image (top left), models (middle), and residuals (right) for Mrk 142. The upper panels display the "simple" models and their residuals, and the bottom 
panels display the "optimal" models and their residuals. The scale of the black bar in the top-left panel is 10 arcsec and the compass in the bottom corner of 
the panel shows the directions north and east. The galaxy image and models are displayed with an inverted logarithmic stretch, and the residuals are displayed 
with an inverted linear stretch centered around zero counts. The bottom-left panel shows the one-dimensional surface brightness of the galaxy (data points), the 
best-fit model (solid line), and each of the individual best-fit model components (PSF = dotted line, all others = dashed lines), with the ellipticity of the galaxy 
displayed below. 



FIG. 3.— Same as Figure 2, but for SBS 1 1 16+583A. 



FIG. 4. — Same as Figure 2, but for Arp 151. The "O" shapes in the images are reflections in the optics from a nearby bright object. 



FIG. 5.— Same as Figure 2, but for Mrk 1310. 



FIG. 6.— Same as Figure 2, but for Mrk 202. 



FIG. 7.— Same as Figure 2, but for NGC 4253 (Mrk 766). 



FIG. 8.— Same as Figure 2, but for NGC 4748. 



FIG. 9. — Same as Figure 2, but for Mrk 290. The small jump in surface brightness at 9" is due to the bright field star in the image. 



Fig. 10.— 



Same as Figure 2, but for NGC 6814. 
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FIG. 1 1. — Top: H/3 BLR radius versus the 5 100 A AGN luminosity. The solid line is the best fit to the data and the grayscale region shows the range allowed 
by the uncertainties on the best fit. The left panel displays all 7 1 datapoints included in this analysis, where the open circles are the new measurements that we 
include for the first time. The right panel shows the fit with Mrk 142 removed, an adopted lag for PG2130+099 of 31 ±4days, and a reddening correction of 
0.26 dex for NGC3227 (see the text for details). The slope does not change appreciably with these adjustments, but the scatter is significantly reduced from 
0. 19 dex to 0. 13 dex. All measurements are plotted with their associated uncertainties, but the error bars are sometimes smaller than the plot symbols. Bottom: 
Residuals of the estimated BLR radii compared to the measured BLR radii using the best fit to the Rblr-L relationship. The dotted lines are Gaussian functions 
with a width equal to the variance in the scatter determined from the best fit, demonstrating the relative normality of the residual distribution. 
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a Distances were estimated from the redshifts of the AGNs except for five cases: NGC 3227, NGC 3783, NGC 4051, 
NGC4151, and NGC 4593. The distances and uncertainties for the last four of these objects are from an average of 
distance moduli to neighboring galaxies in the Tulty et al. (2008) study of the local velocity anomaly, with distance 
uncertainties estimated at 20%. The distance to NGC 3227 is based on the SBF method ITonrv et alJ2001l) for NGC 3226, 
with which NGC 3227 is currently interacting. Additional details are given in the text. 
b Values are from the Schlafiy & Finkbeiner 1201 1) recalibration of the Schlesel et al. i 1998) dust map. 



TABLE 2 
HST Observation Log 
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Surface Brightness Decomposition — Mrk 142 
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NOTE. — Values in square brackets were held fixed during the surface brightness model fitting. 
The STmag magnitude system is based on the absolute physical flux per unit wavelength. 



TABLE 4 

Surface Brightness Decomposition — SBS 1 1 16+583A 
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NOTE. — Values in square brackets were held fixed during the surface brightness model fitting. 
a The STmag magnitude system is based on the absolute physical flux per unit wavelength. 
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TABLE 5 

Surface Brightness Decomposition — Arp 1 5 1 
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10 


sersic 


-15.61 


-10.93 


18.26 


0.50 


1.9 


0.87 


-55.0 






11 


sersic 


-23.94 


-16.20 


17.21 


11.79 


[1.0] 


0.70 


78.3 






12 


sersic 


-31.16 


-19.91 


17.31 


7.23 


[1.0] 


0.36 


-61.5 






13 


PSF 


-33.40 


-20.88 


23.97 














14 


sersic 


-33.43 


-20.84 


18.65 


9.73 


4.3 


0.32 


26.2 






15 


sersic 


-44.52 


-2.88 


19.77 


1.72 


[1.0] 


0.27 


14.8 






merit 




X 2 = 22599108.0 




AW = 10751120 


AW=75 


Xl =2.102 





NOTE. — Values in square brackets were held fixed during the surface brightness model fitting. 
a The STmag magnitude system is based on the absolute physical flux per unit wavelength. 



Bentz, et al. 



TABLE 6 

Surface Brightness Decomposition — Mrk 1310 



Fit 


# 


PSF+sky 


A* (") 


Ay (") 


W^stmag' 




Sky (cts) 


(10 4 cts) 


^ (KT 4 cts) 


Note 






sersic 


Ax (") 


Ay (") 




r e (") 


n 


b/a 


PA (deg) 








sersic3 


Ax(") 


Ay(") 


s b 


n>(") 


n 


b/a 


PA (deg) 








radial 


A* (") 


Ay(") 




1>reak (") 


Ar soft (") 


b/a 


PA (deg) 








power 




r-m (") 


rout (") 


6l ro t (deg) 


a 


^incl (deg) 


8 sky (deg) 








fourier 




mode: a m , 


4> (deg) 


mode: a m 


, <f> (deg) 


mode: a m 


, 4> (deg) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0) 


(8) 


(9) 


(10) 


(11) 


Optimal 


1,2 
3 
4 


PSF+sky 

PSF 

PSF 




-0.03 
-0.15 



0.02 
0.00 


17.42 
19.45 
21.70 




[73.5] 


5.9 


-3.2 






5 


sersic 


-0.02 


0.00 


16.15 


3.29 


4.8 


0.80 


163.0 








fourier 




1: -0.075 


70.9 


3: -0.0086 


-42.9 


4: 0.0092 


6.5 








fourier 




5: -0.0050 


21.7 


6: 0.0093 


-26.8 










6 


sersic 


-0.053 


0.014 


15.34 


4.84 


[1-0] 


0.86 


-99.9 








power 




1.73 


J.J / 


365.9 


0.050 


Jz. I 


"70 O 








fourier 




1: 0.0063 


-40. J 














7 


sersic3 


1.69 


-1.83 


19.41 


7.41 


[1.0] 


0.25 


38.2 








radial,inner 


-0.33 


-45.12 




1.57 


2.76 


0.52 


-37.6 








fourier 




1: 0.56 


69.4 


3:-0.07 


-57.3 


4: 0.093 


14.2 








fourier 




5: -0.062 


-3.7 
















radial,outer 


0.03 


-47.98 




1.77 


0.99 


0.97 


38.3 








fourier 




1: -0.66 


-49.3 


3:0.02 


45.4 


4: 0.073 


34.7 








fourier 




5: -0.047 


18.4 














8 


sersic 


3.50 


-7.55 


22.21 


0.66 


1.3 


0.26 


-31.0 


small galaxy 






fourier 




1:-0.14 


117.2 












merit 




X 2 = 63558184.0 




AW = 12919111 


AW=85 


Xl =4.920 




Simple 


1,2 


PSF+sky 








17.30 




72.5 


6.1 


-5.1 






3 


sersic 


-0.03 


-0.01 


17.93 


0.32 


2.4 


0.89 


-18.7 


bulge 




4 


sersic 


0.06 


0.03 


15.37 


3.39 


0.9 


0.77 


-20.2 


bulge+ring 




5 


sersic 


-0.21 


-0.67 


15.81 


12.03 


[1.0] 


0.71 


-7.2 


disk 




merit 




X 2 = 63561608.0 




AW = 12919172 


AW=19 


Xl =4.920 





NOTE. — Values in square brackets were held fixed during the surface brightness model fitting. 
The STmag magnitude system is based on the absolute physical flux per unit wavelength. 



Low-Luminosity End of the R-L Relationship 



TABLE 7 

Surface Brightness Decomposition — Mrk 202 



Fit 


# 


PSF+sky 


Ax (") 


Ay (") 


^stmag 




Sky (cts) 


ST ( l0 cts > 


d Skv / 1 r\— A . \ 

(10^ cts) 


Note 






sersic 


Ax (") 


Ay (") 


^stmag 


r e (") 


n 


b/a 


PA (deg) 








sersic3 


Ax (") 


Ay(") 


s b 


n,(") 


n 


b/a 


PA (deg) 








radial 


Ax (") 


Ay (") 




'"break (") 


Ar soft (") 


b/a 


PA (deg) 








power 




r m (") 


rout (") 


dmi (deg) 


a 


flincl (deg) 


6» sky (deg) 








fourier 




mode: a m . 


4> (deg) 


mode: a m 


, <f> (deg) 


mode: a m 


, 4> (deg) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Optimal 


1,2 


PSF+sky 








18.13 




34.4 


5.8 


-7.8 






3 


sersic 


0.01 


0.02 


16.47 


0.48 


2.2 


0.79 


-80.3 








fourier 




1: 0.11 


74.3 


3: 0.019 


-6.2 


4: 0.0054 


14. 1 








fourier 




5: 0.0074 


0.9 


6: -0.0064 


6.5 










4 


sersic3 


1.25 


2.74 


9.34 


0.61 


0.8 


0.61 


-28.1 








radial,inner 


-0.60 


-46.35 




1.65 


0.80 


0.66 


-99.1 








fourier 




1: 0.51 


159.9 


3: -0.17 


40.5 


4: 0.13 


27.3 








fourier 




5: 0.074 


1 A 1 

-14. J 


6: -0.041 


-16.5 












radial,outer 


-3.03 


-46.57 




5.06 


14.97 


0.25 


-54.5 








fourier 




1: 9.89 


10.8 


3: -0.39 


24.5 


4: 0.047 


29.6 








fourier 




5: 0.39 


-30.9 


6: 0.54 


26.3 










5 


sersic 


-2.47 


-0.48 


18.48 


1.30 


0.3 


0.72 


45.2 








fourier 




1: 0.67 


35.5 


3: 0.16 


-56.2 


4: O.ll 


35.4 








fourier 




5: 0.098 


20.3 


6: 0.037 


16.8 










6 


sersic 


-0.30 


0.41 


15.66 


16.24 


[1.0] 


0.40 


20.6 








power 




1.18 


2.38 


-206.8 


0.0007 


65.7 


-35.6 








fourier 




1: 0.079 


-96.5 


3: -0.017 


13.1 


4: -0.0058 


7.3 








fourier 




5: 0.0030 


-11.5 


6: 0.0035 


-13.4 










merit 




X 2 = 86674384.0 




N ioi = 12917981 


AW= I0l 


Xl= 6.710 




Simple 


1,2 


PSF+sky 








18.03 




34.3 


5.9 


-8.3 






3 


sersic 


0.02 


0.05 


17.93 


0.23 


1.2 


0.66 


-84.6 


bulge 




4 


sersic 


-0.01 


-0.10 


15.14 


4.79 


4.4 


0.76 


-85.7 


bulge+ring 




5 


sersic 


4.89 


4.86 


18.95 


7.54 


[1.0] 


0.30 


62.72 


disk 




merit 




X 2 = 87365408.0 




AW= 12918056 


AW= 26 


xl = 6.763 





NOTE. — Values in square brackets were held fixed during the surface brightness model fitting. 
The STmag magnitude system is based on the absolute physical flux per unit wavelength. 



Bentz, et al. 



TABLE 8 

Surface Brightness Decomposition — NGC4253 



Fit 


# 


PSF+sky 


Ax (") 


Ay (") 


'^stmag ' 




Sky (cts) 


(ur 4 cts) 


2^ (10 4 cts) 


Note 






sersic 


Ax (") 


Ay(") 


^stmag 


r, (") 


n 


b/a 


PA (deg) 








sersic3 


A* (") 


Ay(") 




rb(") 


n 


b/a 


PA (deg) 








power 




r m (") 


r out (") 


d mt (deg) 


a 


^ind (deg) 


#sky (deg) 








radial 


A* (") 


Ay (") 




r break (") 


Ar soft (") 


b/a 


PA (deg) 








fourier 




mode: a m 


■ 4> (deg) 


mode: a m 


, 4> (deg) 


mode: a m 


, 4> (deg) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Optimal 


1,2 


PSF+sky 








15.69 




28.5 


2.8 


-5.4 






3 


sersic 


-0.12 


0.01 


17.01 


0.22 


3.9 


0.77 


-4.9 








fourier 




1: -0.25 


-36.1 


3: 0.25 


8.2 


4: 0.091 


-30.1 






4 


sersic 


-3.73 


-3.76 


16.82 


3.71 


0.3 


0.16 


-44.7 








fourier 




1: 0.34 


-139.3 


3: -0.43 


3.9 


4. 0.0/6 


43.2 








fourier 




5: 0.073 


O A 
5.4 


6: 0.13 


4.6 










5 


sersic 


-0.04 


0.33 


15.38 


6.47 


[1.0] 


0.48 


30.4 








power 




[0.00] 


12.75 


199.9 


-1.27 


58.4 


50.0 








fourier 




1: 0.26 


140.8 














6 


sersic 


-0.01 


0.10 


13.93 


9.88 


[1.0] 


o cn 
0.39 


An o 
-4/. 9 






7 


sersic3 


2.52 


-14.53 


21.79 


[20.0] 


[1.0] 


0.72 


11.6 








radiafinner 


-0.05 


-45.17 




22.20 


17.44 


0.89 


78.2 








fourier 




1: 0.11 


0.0 


3: -0.14 


-6.7 


4: 0.074 


-28.8 








radial, outer 


6.24 


-39.60 




0.00 


62.86 


0.39 


78.9 








fourier 




1: 0.36 


-130.4 


3: 0.13 


-8.4 


4: -0.11 


18.1 






merit 




X 2 = 18675714.0 




N iot = 10749147 


AW= 84 


xl = 1-737 




Simple 


1,2 


PSF+sky 








15.72 




27.6 


3.9 


-10.8 






3 


sersic 


-0.04 


0.04 


17.62 


0.19 


0.014 


0.62 


-63.0 


nucleus 




4 


sersic 


-0.07 


0.16 


16.46 


1.42 


1.1 


0.56 


-34.5 


bulge 




5 


sersic 


-0.24 


-0.19 


14.68 


7.82 


[1.0] 


0.30 


-48.9 


bar 




6 


sersic 


0.78 


-0.92 


13.51 


16.17 


[1.0] 


0.84 


-74.8 


disk 




merit 




X 2 = 21660426.0 




N iof = 10749197 


AW= 32 


Xl =2.015 





NOTE. — Values in square brackets were held fixed during the surface brightness model fitting. 
a The STmag magnitude system is based on the absolute physical flux per unit wavelength. 



Low-Luminosity End of the R-L Relationship 



TABLE 9 

Surface Brightness Decomposition — NGC 4748 



Fit 


# 


PSF+sky 


Ax (") 


Ay (") 


Wstmag 




Sky (cts) 


2iS (10^ cts) 


(io 4 cts) 


Note 






sersic 


Ax (") 


Ay (") 


^stmag 


r e (") 


n 


b/a 


PA (deg) 








sersic3 


Ax (") 


Ay (") 


Sb 


r b (") 


n 


b/a 


PA (deg) 








power 




fin (") 


fW (") 


9 m (deg) 


a 


9 inc ] (deg) 


9 Ay (deg) 








radial 


Ax (") 


Ay (") 




'"break (") 


Ar sof , (") 


b/a 


PA (deg) 








fourier 




mode: a m 


, <j> (deg) 


mode: a m 


, cf> (deg) 


mode: a m 


, 4> (deg) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(J) 


(8) 


(9) 


(10) 


(11) 


Optimal 


1,2 


PSF+sky 








16.58 




[38.5] 


0.0 


-0.1 






3 


sersic 


-0.02 


0.01 


16.17 


0.19 


4.8 


0.83 


-4.7 








fourier 




1: -0.011 


-115.2 


3: 0.11 


45.3 


4: -0.024 


2.9 








fourier 




5: 0.017 


1.4 


6: -0.008 


7.3 










4 


sersic3 


-0.06 


0.43 


20.96 


[0.80] 


[1.0] 


0.61 


-87.1 








radial,inner 


0.23 


-46.05 




2.62 


2.90 


0.42 


-83.6 








fourier 




1: 0.55 


43.1 


3: -0.14 


-17.0 


4: 0.063 


30.4 








fourier 




5: 0.038 


-25.4 


6: 0.059 


-14.8 












radial,outer 


0.66 


-47.02 




1.21 


0.38 


0.44 


-30.9 








fourier 




1: 0.53 


-73.6 


3: 0.18 


24.3 


4: 0.072 


-22.4 








fourier 




5: 0.015 


-10.3 


6: 0.034 


-27.7 










5 


sersic 


0.26 


-0.18 


16.28 


2.35 


0.5 


0.93 


-75.7 








fourier 




1: 0.14 


15.8 


3: 0.021 


52.2 


4: 0.016 


11.6 








fourier 




5: 0.010 


21.3 


6: -0.0050 


-0.7 










6 


sersic 


0.78 


0.06 


13.86 


17.53 


[1.0] 


0.58 


17.6 








power 




16.76 


18.59 


-22.0 


0.7 


59.3 


10.0 








fourier 




1: -0.13 


75.0 


3: 0.064 


10.6 


4: 0.038 


-21.0 








fourier 




5: 0.026 


22.5 














7 


sersic 


11.79 


-10.61 


17.51 


1.12 


0.7 


0.78 


-83.7 


large 






fourier 




1: -0.28 


A 1 C 


3: 0.029 


20.5 


4. -U.U4U 


-17.7 


companion 




5 


sersic 


12.05 


-10.26 


1 J.OJ 


7.70 


6.2 


U.5Z 


en 1 
-JV. 1 






fourier 




1: -0.093 


29.4 


3: 0.036 


53.2 


4: 0.030 


-12.6 






9 


sersic 


14.05 


-17.73 


15.98 


7.53 


[1.0] 


1.00 


21.1 








fourier 




1: -0.58 


-71.6 


3: 0.19 


10.7 


4: 0.10 


17.4 






10 


sersic 


12.26 


-10.92 


15.36 


40.25 


[1.0] 


0.09 


4.6 








power 




45.40 


61.01 


-5.1 





80.6 


177.9 








fourier 




1: 0.37 


28.2 


3: 0.036 


-18.0 


4: -0.034 


6.2 








fourier 




5: 0.021 


3.9 














1 1 


sersic 


17.47 


-18.36 


19.74 


0.39 


2.9 


0.63 


18.8 


small 






fourier 




1: -0.032 


-47.8 


3: 0.025 


-4.7 


4: -0.012 


18.3 


companion 




merit 




X 2 = 93724888.0 




N io{ = 10749060 


AW= 167 


Xl = 8.719 




Simple 


1,2 


PSF+sky 








16.75 




38.0 


-1.7 


-2.9 






3 


sersic 


-0.02 


0.01 


16.14 


0.14 


8.3 


0.85 


-16.2 


nucleus 




4 


sersic 


0.36 


0.34 


17.36 


0.70 


0.1 


0.78 


74.2 


lens 




5 


sersic 


0.06 


-0.20 


14.70 


5.82 


2.3 


0.76 


73.1 


bulge 




6 


sersic 


-0.20 


0.021 


13.99 


16.80 


[1.0] 


0.69 


-78.7 


disk 




7 


sersic 


12.05 


-10.25 


18.04 


0.50 


3.0 


0.76 


-48.5 






8 


sersic 


12.09 


-10.37 


16.93 


1.16 


0.7 


0.85 


-81.8 






9 


sersic 


11.61 


-9.56 


14.97 


7.66 


[1.0] 


0.64 


23.7 






10 


sersic 


17.47 


-18.36 


19.90 


0.31 


2.4 


0.64 


18.4 






merit 




X 2 = 95017600.0 




Af dof = 10749170 


Afe cc =51 


Xl = 8-840 





NOTE. — Values in square brackets were held fixed during the surface brightness model fitting. 
a The STmag magnitude system is based on the absolute physical flux per unit wavelength. 



Bentz, et al. 



TABLE 10 

Surface Brightness Decomposition — Mrk 290 



Fit 


# 


PSF+sky 


A*(") 


Ay(") 


Wstmag 3 




Sky (cts) 


^ (10- 4 cts) 


^ (lO" 4 cts) 


Note 






sersic 


A* (") 


Ay (") 


si mag 


(") 


n 


b/a 


PA (deg) 








fourier 




mode: a m , 


<p (deg) 


mode: a m 


, <j> (deg) 


mode: a m 


, 4> (deg) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Optimal 


1 

1,^ 
3 
4 
5 


ror+sKy 
PSF 
PSF 
PSF 




-0.16 
0.02 
-0.14 




-0.00 
0.03 
-0.11 


1 J. J4 

19.39 
16.76 
19.77 




37.2 


n i 
i.j 


o n 

-y.u 






6 


PSF 


-0.25 


8.38 


17.79 










star 




7 


sersic 


0.04 


-0.15 


16.48 


1.25 


2.3 


0.85 


-78.2 








fourier 




1: -0.21 


90.3 


3: 0.013 


-25.6 


4: 0.029 


28.6 








fourier 




5: 0.0093 


-5.9 














8 


sersic 


0.083 


-1.13 


15.73 


6.35 


[1.0] 


0.83 


-84.7 








fourier 




1: -0.11 


90.7 


3: -0.010 


7.9 


4: 0.0072 


-19.7 








fourier 




5: 0.013 


4.4 














merit 




X 2 = 28134014.0 




AW = 12856402 


AW=47 


xl = 2.188 




Simple 


1,2 


PSF+sky 








15.32 




37.1 


7.4 


-9.7 






3 


sersic 


-0.04 


0.02 


17.71 


0.16 


0.0 


0.06 


-58.4 


nucleus 




4 


sersic 


-0.03 


-0.07 


15.65 


3.33 


4.3 


0.92 


-79.2 


bulge 




5 


sersic 


0.07 


-0.07 


16.32 


6.92 


[1.0] 


0.81 


-84.9 


disk 




6 


PSF 


-0.25 


8.38 


17.79 










star 




merit 




X 2 = 29751460.0 




AW = 12856420 


AW=29 


Xl= 2-314 





NOTE. — Values in square brackets were held fixed during the surface brightness model fitting. 
a The STmag magnitude system is based on the absolute physical flux per unit wavelength. 



TABLE 11 

Surface Brightness Decomposition — NGC 6814 



Fit 


# 


PSF+sky 


Ax (") 


Ay(") 


W2stmag a 




Sky (cts) 


^ (10-* cts) 


^(HH cts) 


Note 






sersic 


Ax (") 


Ay (") 


^stmag 


r e (") 


n 


b/a 


PA (deg) 








power 




In (") 


roul CO 


0rot (deg) 


a 


0inci (deg) 


sky (deg) 








fourier 




mode: a m . 


. 4> (deg) 


mode: a m 


, 4> (deg) 


mode: a m 


■ <t> (deg) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Optimal 


1,2 
3 
4 


PSF+sky 

PSF 

PSF 



0.01 
-0.07 



0.06 
0.03 


16.52 
17.99 
18.90 




49.4 


14.7 


19.8 






5 


sersic 


-0.17 


0.01 


15.87 


0.90 


1.5 


0.71 


172.6 








fourier 




1: 0.40 


-74.3 


3: 0.090 


40.3 


4: 0.049 


15.7 






6 


sersic 


-0.98 


0.10 


14.34 


4.10 


1.2 


0.83 


67.7 








fourier 




1: 0.21 


-2.8 


3: 0.018 


1.2 


4: 0.026 


-20.0 






7 


sersic 


-7.57 


-0.10 


12.80 


28.42 


[1.0] 


0.66 


-151.3 








power 




-47.98 


148.17 


-386.2 


-1.038 


25.1 


-73.0 








fourier 




1: 0.44 


15.2 


3: 0.15 


-9.4 


4: 0.18 


1.7 






8 


sersic 


4.24 


0.07 


11.72 


44.75 


[1.0] 


0.71 


-39.1 








power 




2.55 


55.83 


-341.7 


0.536 


().() 


-14.5 








fourier 




1: 0.31 


-102.6 


3: 0.083 


49.0 


4: -0.094 


-18.6 






merit 




X 2 = 15661589504.0 




iV dof = 16832317 


AW=74 


X 2 = 930.448 




Simple 


1,2 
3 
4 


PSF+sky 

PSF 

PSF 



0.01 
-0.08 



0.06 
0.03 


16.53 
17.99 
18.91 




45.7 


7.2 


-3.7 






5 


sersic 


0.49 


0.47 


17.31 


1.59 


1.3 


0.49 


73.2 


bulge 




6 


sersic 


-0.08 


-0.01 


15.06 


1.72 


2.1 


0.97 


21.3 


bulge 




7 


sersic 


-0.03 


-0.41 


14.76 


5.87 


0.6 


0.63 


83.6 


bar 




8 


sersic 


-1.21 


0.98 


11.21 


44.47 


[1.0] 


0.98 


84.7 


disk 




merit 




X 2 = 15663450112.0 




AW = 16832352 


AW=39 


X 2 =930.556 





NOTE. — Values in square brackets were held fixed during the surface brightness model fitting. 
d The STmag magnitude system is based on the absolute physical flux per unit wavelength. 
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TABLE 12 

Ground-Based Monitoring Apertures and Measured Flux Densities 



Object Ref. a / obs [5100 A(l +z)] Aperture PA photflam /5100//HST / ga i[5100A(l+z)] 











( 


II x 


) 


(°) 


10 1 ^ pmc cm ^ A p 1 

L \J Clgd Llll T\. C 








Mrk335 


1 


7.683±0.151 


5 





X 


7.6 


90.0 


5.838 





847 


1.559 ±0.078 




1 


8.809 ±0.176 


5 





X 


7.6 


90.0 


5.838 





847 


1.559 ±0.078 




2 


7.490 ±0.268 


5 





X 


12.0 


0.0 


5.838 





847 


1.646 ±0.082 


PG 0026+129 


3 


2.690 ±0.060 


10 





X 


13.0 


42.0 


76.912 





995 


0.379 ±0.019 


PG 0052+251 


3 


2.070 ±0.065 


10 





X 


13.0 


153.4 


5.839 





976 


0.682 ±0.034 


Fairall 9 


4 


5.950±0.085 


4 





X 


9.0 


0.0 


5.791 





888 


2.997 ±0.150 


Mrk59() 


1 


7.895 ±0.170 


5 





X 


7.6 


90.0 


5.796 





848 


3.965 ±0.198 




1 


5.331 ±0.124 


5 





X 


7.6 


90.0 


5.796 





848 


3.965 ±0.198 




1 


6.366±0.137 


5 





X 


7.6 


90.0 


5.796 





848 


3.965 ±0.198 




1 


8.429 ±0.200 


5 





X 


7.6 


90.0 


5.796 





848 


3.965 ±0.198 


3C120 


1 


4.300 ±0.108 


5 





X 


7.6 


90.0 


5.796 





819 


0.624 ±0.031 






3.370 ±0.084 


5 





X 


12.0 


0.0 


5.796 





819 


0.663 ±0.033 


Aknl20 


1 


10.365 ±0.220 


5 





X 


7.6 


90.0 


5.841 





846 


5.549 ±0.277 




1 


7.823 ±0.148 


5 





X 


7.6 


90.0 


5.841 





846 


5.679 ±0.284 


Mrk79 


1 


6.957±0.154 


5 





X 


7.6 


90.0 


5.842 





831 


1.421 ±0.071 




1 


8.487 ±0.156 


5 





X 


7.6 


90.0 


5.842 





831 


1.421 ±0.071 




1 


7.402 ±0.164 


5 





X 


7.6 


90.0 


5.842 





831 


1.421 ±0.071 


PG 0804+761 




5.480 ±0.073 


10 





X 


13.0 


315.6 


5.840 





965 


0.664 ±0.033 


MrkllO 


1 


3.454 ±0.074 


5 





X 


7.6 


90.0 


5.803 





869 


0.665 ±0.033 




1 


3.964±0.081 


5 





X 


7.6 


90.0 


5.803 





869 


0.665 ±0.033 




1 


2.639 ±0.078 


5 





X 


7.6 


90.0 


5.803 





869 


0.665 ±0.033 


PG 0953+414 


3 


1.560 ±0.032 


10 





X 


13.0 


31.7 


5.841 


1 


099 


0.224 ±0.011 


NGC 3227 


5 


11.617±0.109 


5 





X 


7.5 a 


0.0 


5.800 





814 


6.838 ±0.631 


Mrk 142 


6 


2.050 ±0.045 


4 





X 


9.4 a 


90.0 


4.632 





916 


0.650 ±0.032 


NGC 3516 


5 


20.957 ±0.248 


5 





X 


12.0 a 


0.0 


5.828 





816 


15.36± 1.54 


SBS 1116+583A 


6 


1.088 ±0.048 


4 





X 


9.4 a 


90.0 


4.632 





899 


0.957 ±0.048 


Arpl51 


6 


1.835 ±0.079 


4 





X 


9.4 a 


90.0 


4.632 





000 


1.241 ±0.062 


NGC 3783 


7 


11.380 ±0.604 


5 





X 


10.0 


0.0 


5.795 





801 


4.717 ±0.472 


Mrk 1310 


6 


1.870 ±0.073 


4 





X 


9.4 a 


90.0 


4.632 





000 


1.502 ±0.075 


NGC 4051 


5 


13.667 ±0.194 


5 





X 


12.0 a 


0.0 


5.799 





814 


8.738 ±0.874 


NGC4151 


8 


23.800 ±0.539 


5 





X 


12.0 a 


o.o b 


5.801 





811 


17.03 ± 1.70 


Mrk 202 


6 


1.698 ±0.053 


4 





X 


9.4 a 


0.0 


4.632 





889 


1.395 ±0.070 


NGC 4253 


6 


4.590 ±0.092 


4 





X 


9.4 a 


60.0 


4.632 





877 


2.920 ±0.146 


PG 1226+023 


3 


21.300 ±0.579 


10 





X 


13.0 


171.2 


5.845 





971 


1.303 ±0.065 


PG 1229+204 


3 


2.150±0.039 


10 





X 


13.0 


291.5 


5.842 





915 


1.326 ±0.066 


NGC 4593 


9 


15.853 ±0.319 


5 





X 


12.75 


0.0 b 


5.830 





820 


7.837 ±0.743 


NGC 4748 


6 


4.360 ±0.075 


4 





X 


9.4 a 


0.0 


4.632 





874 


3.207 ±0.160 


PG 1307+085 


3 


1.790 ±0.040 


10 





X 


13.0 


186.5 


76.912 





990 


0.230 ±0.012 


IC4329A 


10 


5.790 ±0.297 


5 





X 


10.0 


90.0 


5.831 





824 


3.570 ±0.357 


Mrk 279 


11 


6.897 ±0.492 


5 





X 


7.5 


90.0 


5.796 





861 


2.918±0.146 


PG 141 1+442 


3 


3.710±0.054 


10 





X 


13.0 


347.0 


5.842 





950 


0.769 ±0.038 


NGC 5548 


12 


9.916 ±0.494 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


7.252 ±0.353 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


9.396 ±0.380 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


6.720 ±0.296 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


9.062 ±0.351 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


9.760 ±0.424 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


12.091 ±0.312 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


10.563 ±0.407 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


8. 120 ±0.306 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


13.468 ±0.548 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


11.832 ±0.444 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


6.981 ±0.304 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




12 


7.032 ±0.334 


5 





X 


7.5 


90.0 


5.801 





833 


3.752 ±0.375 




13 


6.630 ±0.243 


5 





X 


12.75 


0.0 b 


5.801 





833 


4.341 ±0.434 




6 


6.121 ±0.130 


4 





X 


9.4 a 


60.0 


5.801 





833 


3.537 ±0.354 




5 


6.766 ±0.086 


5 





X 


12.0 a 


0.0 


5.801 





833 


4.273 ±0.427 


PG 1426+015 


3 


4.620 ±0.067 


10 





X 


13.0 


341.4 


76.912 





957 


1.134±0.105 


Mrk 817 


1 


6.098 ±0.120 


5 





X 


7.6 


90.0 


5.796 





864 


1.489 ±0.142 




1 


4.998 ±0.104 


5 





X 


7.6 


90.0 


5.796 





864 


1.489 ±0.074 




1 


5.008 ±0.108 


5 





X 


7.6 


90.0 


5.796 





864 


1.489 ±0.074 




5 


6.896 ±0.070 


5 





X 


12.0 a 


0.0 


5.796 





864 


1.723 ±0.086 


Mrk 290 


5 


2.575 ±0.035 


5 





X 


12.0 a 


0.0 


4.632 





898 


1.342 ±0.067 


PG 1613+658 


3 


3.490 ±0.044 


10 





X 


13.0 


164.2 


5.842 





977 


1.460 ±0.073 


PG 1617+175 


3 


1.440 ±0.021 


10 





X 


13.0 


253.0 


76.912 





985 


0.336 ±0.017 


PG 1700+518 


3 


2.200 ±0.011 


10 





X 


13.0 


183.5 


5.842 


1 


403 


0.338 ±0.017 


3C 390.3 


14 


1.732±0.188 


5 





X 


7.5 


90.0 


5.801 





899 


0.826 ±0.041 




15 


6.630 ±0.074 


3 





X 


9.4 


90.0 


5.801 





899 


0.725 ±0.036 


NGC 6814 


6 


6.470 ±0.238 


4 





X 


9.4 a 


150.0 


4.632 





834 


4.132 ±0.207 


Mrk 509 


1 


10.920 ±0.231 


5 





X 


7.6 


90.0 


76.912 





887 


2.434 ±0.122 


PG 2130+099 


2 


3. 100 ±0.076 


5 





X 


12.0 


0.0 


5.794 





913 


0.574 ±0.029 


NGC 7469 


16 


19.023 ±0.572 


5 





X 


12.0 


0.0 


5.836 





822 


8.219 ±0.822 



REFERENCES. —1. Peterson et al. i 1998). 2. Grier et al. 12012a), 3. Kaspi et al. (2000), 4. Santos-Lleo et al. 1 1997), 5. Denney et al. 1 2009), 6. Bentz et al. 
(2009b), 7. Stirpe et al. 1 1994), 8. Bentz et al. 1 2006b), 9. Denney et al. 1 2006), 10. Winge et al. 1 1996), 11. Santos-Lleo et al. i 2001 ), 12. Peterson et al. 1 2002) 
and references therein, 13. IBentz et al.H2007l) , 14. IDietrich et alj 119981) ' 15. IDietrich et alj|2012b , 16. |peterson et alj|2013b . 

NOTE. — Flux densities are tabulated at rest-frame 5 100 A in fx units of ergs s" 1 cm" 2 A -1 . Position angle is defined as degrees east of north. Photflam is the 
inverse sensitivity of the HST detector, listed here in units of 10~ 19 ergs cm" 2 A -1 electron -1 . 
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TABLE 13 

Rest-frame H/3 Time Lags and 5 1 00 A Luminosities 



Object 



H/3 Time Lag 
(days) 



/agn [5100A(1+z)] AL a , agn (5100A) 



/in-15 -1 -2 a-1\ 
(10 ergs s cm A ) 


/ 1 r\44 — 1 \ 

(10 ergs s ') 


6.12±0.19 


43.70 ±0.06 


7.25 ±0.21 


43.78 ±0.05 


5.84 ±0.29 


43.68 ±0.06 


2.31 ±0.07 


44.91 ±0.02 


1.39 ±0.08 


44.75 ±0.03 


2.95 ±0.23 


43.92 ±0.05 


3.93 ±0.33 


43.53 ±0.07 


1.37 ±0.31 


43.07 ±0.11 


2.40 ±0.31 


43.32 ±0.08 


4.46 ±0.34 


43.59 ±0.06 


3.68 ±0.12 


44.01 ±0.05 


2.71 ±0.10 


43.87 ±0.05 


4.82 ±0.45 


43.92 ±0.06 


2. 14 ±0.43 


43.57 ±0.10 


5.54±0.18 


43.57 ±0.07 


7.07 ±0.19 


43.67 ±0.07 


5.98 ±0.19 


43.60 ±0.07 


4.82 ±0.09 


44.85 ±0.02 


2.79 ±0.09 


43.62 ±0.04 


3.30 ±0.09 


43.69 ±0.04 


1.97 ±0.09 


43.47 ±0.05 


1.34 ±0.04 


45.13 ±0.01 


4.78 ±0.73 


42.24 ±0.11 


1.40 ±0.07 


43.54 ±0.04 


5.60 ±1.74 


42.73 ±0.21 


0.13 ±0.08 


42.07 ±0.28 


0.59 ±0.12 


42.48 ±0.11 


6.66 ±0.80 


42.55 ±0.18 


0.37 ±0.1 3 


42.23 ±0.17 


4.93 ±1.00 


41.96±0.20 


6.77 ±1.98 


42.09 ±0.22 


0.30 ±0.11 


42.20 ±0.1 8 


1.67 ±0.23 


42.51 ±0.13 


20.00 ±0.59 


45.90 ±0.02 


0.82 ±0.10 


43.64 ±0.06 


8.02 ±0.90 


42.87 ±0.1 8 


1.15 ±0.24 


42.49 ±0.13 


1.56 ±0.04 


44.79 ±0.02 


3.98 ±0.53 


43.64 ±0.08 


2.94 ±0.08 


44.50 ±0.02 


6. 16 ±0.65 


43.33 ±0.10 


3.50 ±0.55 


43.08 ±0.11 


5.64 ±0.57 


43.29 ±0.10 


2.97 ±0.51 


43.01 ±0.11 


5.31 ±0.55 


43.26 ±0.10 


6.01 ±0.60 


43.32 ±0.10 


8.34 ±0.52 


43.46 ±0.09 


6.81 ±0.59 


43.37 ±0.09 


4.37 ±0.52 


43.18±0.10 


9.72 ±0.69 


43.52 ±0.09 


8.08 ±0.61 


43.44 ±0.09 


3.23 ±0.52 


43.05 ±0.11 


3.28 ±0.54 


43.05 ±0.11 


2.29 ±0.54 


42.90 ±0.13 


2.58 ±0.42 


42.95 ±0.11 


2.49 ±0.49 


42.93 ±0.12 


3.49 ±0.14 


44.57 ±0.02 


4.61 ±0.20 


43.73 ±0.05 


3.51 ±0.15 


43.61 ±0.05 


3.52 ±0.15 


43.61 ±0.05 


5.17±0.14 


43.78 ±0.05 


1.23 ±0.10 


43.11 ±0.06 


2.03 ±0.11 


44.71 ±0.03 


1.10±0.03 


44.33 ±0.02 


1.86 ±0.03 


45.53 ±0.01 


0.91 ±0.20 


43.62 ±0.10 


5.91 ±0.09 


44.43 ±0.03 


2.34 ±0.38 


42.05 ±0.29 



Mrk335 



PG 0026+129 
PG 0052+251 
Fairall 9 
Mrk59() 



3C 120 
Ark 120 
Mrk79 



PG 0804+761 
MrkllO 



PG 0953+414 
NGC3227 
Mrk 142 
NGC3516 
SBS1116+583A 
Arpl51 
NGC 3783 
Mrk 1310 
NGC 4051 
NGC4151 
Mrk 202 
NGC 4253 
PG 1226+032 
PG 1229+204 
NGC 4593 
NGC 4748 
PG 1307+085 
Mrk 279 
PG 1411+442 
NGC 5548 



PG 1426+015 
Mrk 8 17 



Mrk 290 
PG 1613+658 
PG 1617+175 
PG 1700+518 
3C 390.3 

NGC 6814 



16.8+* 
12.5+" 5 
14 3 
111.0+241 

89.8+1 
17.4+3 

7+3 



20.7_ 
14.0+.! 
29.2+ 
28.8+ 
38.1+? 5 
25.9+^ 
47.1+P 
37.1+|! 
9.0+fl 

16.0+|4 
146.9+Jf 



24.3+ 



3 



20.4+l° 3 5 

-3-7 -J + li9 

-10.0 

150.1+^« 

3 75 +0 4 ?6 
J -' J -Q.S2 

2 74+0.73 
Z - ' -0-83 

11.68 

2.31 

3 99^> M 
10:20+|| 

1 C7+0.54 

6.58+ ;„ 

3.05+} || 

6. 16+™ 

306.8^ 

37.8! f 5 j 
a 7r rf> 75 

162 
5 . 5 5 ^ 22 

105 60+ 3 ° " 

16.70+H8 
124.30+j" 
19.70+1-50 

18.60+Hjj 

15.90+l : |g 

1 1 00+1 90 

1 3 00+ 1 60 
1J - -1.40 

21 70+ 2 - 60 
-2.69 

16 40 +t -20 

-119 

17 <S0+ 2 -°0 
1/ - ,u -1.60 

26.50+ 4 -| 
"'Si 

D - ,u -3.70 

14.30+^0 

6 30+ 2: ™ 
-2-39 

4 | 8 +0.86 

12.40+ 2 -g 
95.0+2 i)! ' 
19 00+ 3 " 90 

* -3.7,0 
, c Q +"0 

-3.50 

33 60+ g ^ 8 

14.04+|« 

8.72+ |g 

40.10+J|;0 

71 50+ 29 - 6 
-33 .70 

251.8+ 4 || 
23 60+*-20 
4fi 40+ 3 - 60 
+0 S fi 



6.64 



'-o-?o 



Low-Luminosity End of the R-L Relationship 



Note 



TABLE 14 
H/9 Rblr-L Fits 



All 

Clean 

Clean2 

ExtCorr 

Clean+ExtCorr 

Clean2+ExtCorr 



71 

70 
69 
71 
70 
69 



I CO7+0.031 
^ -0.031 

i 555+0.024 

■ — 0.024 

1 526 + ™ 
1 554+0 024 
I 5^0+0.02? 



5 3 3+ 035 

"■ JJJ -0.033 
542+ - 027 

0.546^1 
0.549™ 



03 6+ 009 



0.546 
0.537 



-o 
0.018+0 

0.018;? 

0.037™ 

0.018; 

0.018 



pns 
006 



NOTE. — The scatter is reported as a 2 , the variance of the measure- 
ments around the best-fit relationship. All: All individual measurements 
are included in the fit. Clean: The potentially erroneous time-lag measure- 
ment for Mrk 142 is excluded from the fit, and the potentially biased time lag 
for PG2130+099 is replaced with r = 31 ± 4days. Clean!: The lags for 
both Mrk 142 and PG 2130+099 are excluded. ExtCorr: The luminosity of 
NGC 3227 is corrected for 0.26 dex of internal extinction. 



